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Summary 


This paper describes an approach which aims at bridging the gap between the traditional 
Reynolds-averaged Navier-Stokes (RANS) approach and the traditional large eddy 
simulation (LES) approach. This approach affords an intermediate resolution of turbulent 
scales relative to those of RANS and LES, and has the characteristics of the very large 
eddy simulation (VLES). In the present work, the very large scales of turbulence are 
directly calculated, and the effects of the unresolved scales are accounted for by an eddy 
viscosity model plus the nonlinear source terms representing the effects of anisotropy and 
rotation. The dependent variables and governing equations are based on a temporal 
filtering with a constant filter width. The contents of both resolved and unresolved scales 
are regulated by the width of the temporal filter. In the present approach, the dependent 
variables and governing equations will naturally evolve from RANS to VLES and further 
towards LES, when the width of the temporal filter decreases from the turbulent integral 
time scale to its fraction and all the way towards the Taylor micro time scale. We call this 
approach the partially-resolved numerical simulation (PRNS). The main features of the 
PRNS are: the filtered governing equations are mathematically grid invariant; the 
subscale model is the one evolved from state-of-the-art models used in the RANS 
approach; and its nonlinear formulation explicitly accounts for the important features 
occurring in the regime of the very large eddy simulation. 

The concept of PRNS approach, the basic equations and the subscale model are described 
in detail. Systematic simulations using the National Combustion Code (NCC) have been 
carried out for fully developed turbulent pipe flows at different Reynolds numbers to 
evaluate the PRNS approach. The calculated results and the influences of the numerical 
settings used in the NCC are assessed. 

Also presented are the sample results of two demonstration cases: non-reacting flow in a 
single injector flame tube and reacting flow in a Lean Direct Injection (LDI) hydrogen 
combustor. These results indicate that the PRNS approach is able to capture important 
large scale turbulent structures and thereby improve the fidelity of numerical simulations 
while keeping the computational costs much lower than that required by the LES. 


I. Introduction 

1. Background of PRNS 

Flows in combustors are three dimensional and intrinsically unsteady due to massive 
separation and strong swirling. These flows contain large scale, coherent structures which 
play a very important role in determining the turbulent mixing of the fuel, air and other 
reactants. This, in turn, will heavily influence the flame structure and flame stability via 
turbulence-chemistry interaction, consequently, also the efficiency, robustness and 
emissions of the combustion. Apparently, the capability of explicitly capturing these 
coherent flow structures in the computation is critical for improving the fidelity of 
simulations of turbulent combustion. 
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Since, by its definition, the state-of-the-art Reynolds averaged Navier-Stokes (RANS) 
approach has been aimed at establishing the long-time averaged statistics of flow 
variables, it can not provide any information on those instantaneous coherent flow 
structures; hence it is very difficult to accurately account for the turbulent mixing and the 
turbulence-chemistry interaction in the RANS calculations. On the other hand, the state- 
of-the art traditional large eddy simulation (LES), which can capture the unsteady 
coherent flow structures in the computation but typically requires grid sizes down to the 
inertial sub-range of the local turbulence, is often too costly for practical engineering 
applications. There is a need for approaches which are capable of capturing the 
dynamically dominant large-scale flow structures in the computation but at a reasonable 
computational cost. 

One type of these approaches is illustrated by the recent development of various hybrid 
RANS/LES methodology (References [1, 2, 3, 4, 5, 6, 7]). The basic strategy is to reduce 
the turbulent eddy viscosity in the fine grid regions to promote LES type of simulation 
and revert to RANS type of simulation in the coarse grid regions. The switch between 
RANS and LES is usually based on a criterion in terms of local grid spacing and local 
flow quantities. Most of them invoke spatially-filtered equations for the resolved scales 
of turbulence, and the sub-grid model explicitly has the local computational mesh size as 
a parameter. In addition, the local grid spacing is typically considered as the local width 
of the spatial filter. There are several troublesome issues associated with the spatially- 
filtered equations and the use of computational mesh size as a parameter in the sub-grid 
model as discussed in Reference [8], However, if we use a temporal filter with a fixed 
filter width to define the large scales of turbulence, together with a sub-scale model 
which does not have a grid spacing parameter in its constitutive equation, then all of these 
issues can be avoided. This leads to our proposed approach termed as the partially 
resolved numerical simulation (PRNS). 


2. Basic Concept of PRNS 

PRNS is aiming at the very large eddy simulation (VLES) of high Reynolds number 
turbulent flows using a relatively coarse grid which is comparable to the one used in the 
RANS approach. PRNS is based on the concept of using the temporal filter to define the 
large or very large scales of turbulence. The basic equations for the resolved scales of 
turbulence are the temporally filtered Navier-Stokes equations. The large scales of 
turbulence are directly calculated by solving these temporally filtered equations, just like 
any other spatially filtered LES approaches. However, the effects of the unresolved scales 
are modeled by a dynamic equation system which is evolved from the state-of-the-art 
RANS models. The contents of both resolved and unresolved turbulence are notionally 
regulated by a “resolution control parameter” (RCP), which is related to the ratio of the 
temporal filter width to the global time scale of the turbulent flow of interest. In practice, 
the value of the RCP affects the effectual separation of the resolved scales from the 
unresolved ones in the computation. The basic equations of PRNS as well as the subscale 
models are thus, in theory, grid independent or grid invariant. The nature of the 
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mathematical solution of the PRNS equations will mainly depend on the magnitude of 
RCP under the imposed boundary and initial conditions. As RCP varies from 0 to 1, the 
nature of the solution should vary from the direct numerical simulation (DNS), LES, 
VLES and toward RANS. In the actual computational simulations, however, the 
appropriate grid resolution must be provided for a given value of RCP to support a 
meaningful turbulent simulation. Overly fine grid will not change the nature of the 
simulation, although numerical effects such as less numerical dissipation are expected. 
This is because the basic equations and subscale models of PRNS do not involve any 
parameters containing the grid size, unlike the traditional LES approaches that are 
directly dependent on the mesh size. 

The basic equations of PRNS and the subscale models are described in Chapter II, 
Sections 1.2 and 1.4. Some remarks about planning the VLES, choosing an appropriate 
value of RCP, and the main features of the dynamic equation system of the subscale 
model are discussed in Chapter II, Section 1.4.5. 


3. The National Combustion Code (NCC) 

The approach of PRNS can be easily implemented into an unsteady computational fluid 
dynamics (CFD) code. In the present effort, we have used the National Combustion Code 
developed at NASA GRC [9, 10, 11, 12, 13, 14, 15], The NCC is an integrated system of 
modules based on unstructured meshes and running on parallel computing platforms. It 
has five major modules: a gaseous flow module solving the three-dimensional, unsteady, 
compressible Navier-Stokes equations; a turbulence module including the nonlinear k-s 
models; a chemistry module solving the species conservation equations, or using the 
Intrinsic Low Dimensional Manifold (ILDM) method to determine the species mass 
fractions; a turbulence-chemistry interaction module employing eddy-breakup model, or 
assumed pdf approach, or solving the transport equation of the probability density 
function (pdf) for species and enthalpy; and a spray module solving the liquid droplet 
transport equations. The interacting multi-phase, physical-chemical processes embodied 
in these major modules are emulated via a hybrid Eulerian-Lagrangian-Monte Carlo 
solution algorithm. Over the years, the NCC has served as a test bed for assessing 
turbulent combustion models and computational technologies in a semi-engineering 
environment. We have been prudently using the NCC within its current limitations for 
technology program support while concurrently improving its fidelity and extending its 
capability. 

The implementation and evaluation of the PRNS mainly involve the gaseous flow module 
and the turbulence module, therefore, a brief description of these two modules is in order. 
The flow module solves unsteady, three-dimensional, compressible Navier-Stokes 
equations. The discretization begins by dividing the computational domain into a large 
number of mesh elements, which can be of mixed type. A central-difference finite- 
volume scheme is used for spatial discretization while a three-time-level backward 
differencing is used for temporal discretization. Therefore, the nominal discretization 
accuracy is second order both in space and in time. These discretized equations are then 
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advanced temporally to the next physical-time step by a so called “dual- time- step” 
approach, in which a four-stage Runge-Kutta scheme is used to iteratively drive the 
solutions of the “pseudo-time” towards convergence with respect to the pseudo time. 
Upon the convergence of this “inner” iteration, the implicit advancement to the next 
physical-time step is achieved. In order to accelerate the convergence, particularly when 
the flow Mach number is small, the pseudo-time term is preconditioned, and the static 
pressure is decomposed into a constant reference pressure and a gauge (or fluctuating) 
pressure. Sometimes, residual smoothing is also applied to smooth the computed 
residuals to enhance the convergence rate. Since central differencing is employed for the 
spatial discretization, artificial dissipation is added to maintain the numerical stability, we 
have been using the Jameson dissipation operator, which is a blend of second-difference 
and fourth-difference terms. Consequently, for a practical computation, the numerical 
setting will then include the truncation errors of the discretization, the level of the 
reference pressure, the convergence criterion for the inner iteration, the use of residual 
smoothing, and the amount of the added artificial dissipation. A sufficiently “clean” 
numerical setting is the prerequisite for a rigorous evaluation of the turbulence models. 
By clean, we mean that the spurious effects originated from the inherent numerical 
artifacts are much smaller than the effects originated from the employed physical models 
on the solutions. In other words, the behaviors of the calculated solutions are controlled 
by the physical models, not by the numerical artifacts, although some limited yet small 
contaminations are always present. 

The turbulence module is based on the k-s models, and has the following options: high 
Reynolds number wall function or low Reynolds number wall integration, linear or non- 
linear turbulent stress-strain relationship. Details of an advanced nonlinear k-e RANS 
model and a recently developed generalized 3-D turbulent wall function can be found in 
references [15] and [17], respectively. The module also contains some LES sub-grid 
models [18, 19]. 


4. Scope of Assessment Effort 

Recently, PRNS is being implemented in the NCC for performing the very large eddy 
simulation of combustor flows. The focus of the current effort is the fundamental 
assessment of the implemented subscale model and the establishment of an appropriate 
numerical setting. Experiences have shown that both the physics based model and the 
appropriate numerical setting are equally important for yielding an accurate numerical 
simulation. In the past, the NCC has successfully produced many RANS simulations for 
various reacting combustor flows. The code is robust and has a set of default numerical 
parameters to speed up the convergence and to stabilize the numerical simulation. Some 
of these parameters/options (e.g., the reference pressure, residual smoothing, added 
numerical dissipations, etc.) are common practice, and the RANS results are often not 
overly sensitive to a range of their values. However, we have observed that these 
parameters/options can significantly affect and, in some cases, even dominate the 
simulation of turbulent fluctuations. Subsequently, we have conducted systematic studies 
to sort out the effects of each of these parameters/options on the calculation of turbulent 
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fluctuations, and, as a result, we have established a sufficiently “clean” numerical setting 
for the PRNS calculation. The simulations presented in this paper are all conducted using 
this numerical setting as the baseline. 

At this stage of the assessment of the PRNS approach, we have selected the fully 
developed turbulent pipe flows as the focus, as they are well known fundamental flows 
with available experimental data [20]. Furthermore, these cases are amenable to the use 
of periodic boundary conditions in the simulations, thereby, free from the complicating 
issues of inflow and out flow boundary conditions. In this report, we will present the 
assessment results obtained for a moderately large Reynolds number of 150,000 and a 
relatively low Reynolds number of 15,000 (based on the centerline axial velocity and the 
pipe diameter). The basic criteria for a successful simulation of large and very large 
scales of turbulence are set as follows: it must reveal a broadband scales of turbulent 
fluctuations, the power spectrum has to show at least two or three orders energy 
variations from large to small scales, and it must produce reasonable or better mean 
profiles than those produced by the RANS simulation, when compared with experimental 
data. 

Finally, we also present the sample results of two demonstration cases to indicate the 
intended applications of the PRNS approach. These two cases are: non-reacting flow in a 
single injector flame tube and the reacting flow in a Lean Direct Injection (LDI) 
hydrogen combustor. 


II. Partially Resolved Numerical Simulation 

In this chapter, we will first describe the basics of PRNS approach for the very large eddy 
simulation, the advantages of its grid-invariant feature and its subscale dynamic equation 
system. Then, we will present the critical evaluations of PRNS using the NCC code for 
fully developed pipe flows at various Reynolds numbers. 

Some interesting phenomena were observed, which seem to reveal some deficiencies 
resided in all of the turbulent simulations that are based solely on the subscale eddy 
viscosity concept. More advanced subscale models may be needed to more 
comprehensively mimic the physical interactions between the resolved and unresolved 
scales, especially when the unresolved scales become larger or the flows are at lower 
Reynolds numbers. 

In addition, sample results of applications to complex non-reacting and reacting flows 
will also be presented to show the intended applications of this approach. 

1. Basics of Partially Resolved Numerical Simulation 

In the PRNS approach, the resolved turbulent scales are defined by a temporal filter. This 
is based on the premise that small scale motions tend to have small time scales. A 
temporal filter with a fixed small width can then capture the large scale turbulence. 
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The PRNS equations are the temporally filtered Navier-Stokes equations, which govern 
the motion of the resolved large scale turbulence. There are unclosed terms in these 
PRNS equations, which represent the effects from the unresolved small scale turbulence, 
and they must be modeled in order to solve the PRNS equations to obtain the numerical 
solution of large scale turbulence. 

It is noted here that the above procedure in PRNS is very similar to that in the traditional 
LES, except that the filter used to define the large scale turbulence is different in these 
two approaches. PRNS uses a temporal filter while LES uses a spatial filter (with its 
width often being identical to the numerical grid size). There are substantial advantages 
in using a temporal filter. For example, PRNS equations are mathematically grid- 
invariant; they are, in theory, valid for any type of numerical grids (structured or 
unstructured, homogeneous or inhomogeneous); PRNS equations have a unified feature 
facilitating the evolution from DNS, LES and towards RANS when the width of the 
temporal filter varies from zero, a small but finite value and towards a large value (e.g., 
the integral time scale). In addition, some important turbulent correlation relationships 
can be unambiguously established allowing an exact comparison between the 
measurements and numerical simulations. 

In order to be consistent with the grid invariance of the PRNS equations, the subscale 
model of PRNS must also be grid-independent; it can be formed by invoking a more 
general constitutive relationship in conjunction with a dynamic equation system 
consisting of the subscale turbulent kinetic energy and dissipation rate. This type of 
models enables PRNS to perform a very large eddy simulation with a relatively coarse 
grid, hence, consuming less computing resources. 

In a preliminary study, we have noticed that the non-linear constitutive relationship for 
the subscale model exhibits significant advantages over the eddy viscosity model. This 
suggests that the interactions between the resolved and the unresolved scales are not just 
through the eddy viscosity, and that the nonlinear terms are needed as additional sources 
to mimic other types of physical interactions. 


1 . 1 Definition of Large-Scale T urbulence 

Using a homogeneous temporal filter G(t-t ), the large-scale turbulent variable (j) and 
its density-weighted variable (/) can be defined as 

4>(t,x i ) = f (j)(t' ^^Git-t^dt , (1) 

J P 

where the integral is over the entire time domain and G satisfies the normalization 
condition: J G(t-t )dt =1. There are many such temporal filters. One of the simplest 
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filters is the top hat filter. For example, the following top hat filter with the width A r is 
used in the current PRNS: 


G(t-i) = \ 1/Ar ’ if \* 

\ 0 , otherwise. 


Using this filter, the left part of Equation (1) will become 

^ t+Aj* / 2 

0(t,x t ) = — [ Mt,x t )di. 

A J 

T t- A t /2 


( 2 ) 


(3) 


Equation (3) reveals an unified feature of (/> and </> , because they will become the exact 
Reynolds-averaged quantity and Favre-averaged quantity when A r — > go . On the other 
hand, they will become the instantaneous turbulent quantity as A r — » 0 . For a finite A r 5 
they represent the quantities of large scale turbulence. 


1.2 Basic Equations for PRNS 

Performing operation (1) on the Navier-Stokes equations, we obtain a set of basic 
equations for the resolved large scale turbulence ( <j> and <j ) ): 


(M) , + (M“,) =- Pj-Tyj + 2/i s g --Sy/u. 

,J V J 


kk 


~ / ^ 

(p e ) t +(ptii e\, = (kTj) , + P s kk-%i + \ 2 P s ij s ij ~TPhk s u 

p t +(pUi). = 0, p = pR f , 


+ 0 5 


(4) 

(5) 

( 6 ) 


where v =(u ij + u j ^j 2. The symbols ( ) t and ( ) . represent the temporal and spatial 


derivatives, respectively, p , u t ,, T , p , e , and Q are the density, velocity, temperature, 
pressure, internal energy per unit mass, and the radiation rate, p and k are the viscosity 
and heat conductivity. R is the universal gas constant, zv and q t are the extra terms that 
are created during the process of temporally filtering Navier-Stokes equations: 


Ty = p (u t Uj - up J ) , q,=p (u,e - up). 


(V) 


They represent the effects from the unresolved subscale turbulence. They are not in a 
closed form, hence must be modeled. We refer them as the unresolved turbulent stresses 
and heat fluxes. 
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Note that there are two types of resolved large-scale turbulent variables, <j) and <f> , 

appeared in the above equations. Where (f> is defined by a temporal filtering and (j) is 
defined by a density-weighted temporal filtering. 


1.3 Time- Averaging Relationships 

Since various filters can be used to define the calculated, resolved turbulence quantities, a 
practical, yet somewhat overlooked, question is how to conduct an apple-to-apple 
comparison between the calculated results and the experimental data. 


Most of the turbulent experimental data are the “mean” values of turbulent variables (e.g., 
velocity, temperature, pressure, etc.); they are either the pure time-averaged values (for 
incompressible flows) or the density-weighted time-averaged values, i.e., the Favre- 
averaged values (for compressible flows): 


j 7-/2 j 7-/2 

M = I M = J Pt*, 


T — »co f 


( 8 ) 


-77 2 


where ^ represents an instantaneous turbulent quantity, it can be u l , u^iij , or u j u J u k , etc., 

T is the entire time domain, is the pure time-averaged (Reynolds-averaging) quantity 

and \<p\ is the density-weighted time-averaged (Favre-averaging) quantity. Such defined 

experimental data are actually based on the assumption that the measured turbulent flows 
are statistically steady, or at least approximately steady. Otherwise, the experimental data 
of “mean” values must be redefined using the ensemble average of many repeated 
realizations of the same experiment, i.e., 


{,/} = lim — f> |r » , ^T- = lim —] ^ 

j { p } [p]M 




(0 


(9) 


where M is the number of realizations and </> <n and ( p<f) f ] are the individual realization. 
The symbol { } represents the ensemble average. For statistically steady turbulent flows, 
the time average and the ensemble average become identical, i.e. 

M-W. do) 

In summary, the experimental “mean” values for the compressible flow are either the 
density-weighted time average [^] or the density-weighted ensemble average { p<f >\ /{/?}, 
expressed by Equations (8) and (9), respectively. 
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Now, in PRNS, the temporally filtered values </>,</) of turbulent variables are directly 
calculated. As a result, this enables an apple-to-apple comparison between the PRNS data 
and experimental data. In Reference [8], relations between the PRNS data and the 
experimental data have been established for the first order, second order and higher order 
correlations. Furthermore, the relations for the Reynolds stresses and the scalar fluxes 
(which are second order moments formed by the zero-mean fluctuating velocities and 
scalars) can also be derived. Here we list these relations under the statistically steady 
condition: 

[u i ] = {pu i )/(p), [e] = (pe)/{p) (11) 

H = (M £ )/(p). [«,.«,.] = ( 12 ) 

Equation (11) indicates that the experimentally measured mean velocity and mean scalar 
(terms on the left hand side of the equation) can be compared directly with the post- 
processed PRNS data (terms on the right hand side). Equation (12) cannot be used for 
direct comparison, because the unclosed terms are involved; however, they are useful for 
establishing the following relations for the Reynolds stresses and the scalar fluxes: 



(13) 

*,=z; +<?,)/(/>), 

(14) 

where, 

R tj s {p (“/ -l u i])( u j-[ u j]))/(p)= \ u i u j ] - h ] I u j ] ’ 

(15) 

R i=(p ( u i - [“, ]) [ e - [ e ])) /{p) = VA - Vi ] [ e ] > 

(16) 

T v = {PWj ) {P) ~ (PR > {P S J )/ (Pf > 

(17) 

T i={pu i e)l{p)-{pu i ){pe)l(p ) 2 . 

(18) 


The experimentally measured Reynolds stresses R tj and the scalar fluxes R are 

determined according to Equations (15) and (16), respectively. Their PRNS counterparts 
consists of two components: the first component, T. or T., can be determined by using 

the directly calculated variables according to Equation (17) or (18); the second 
component, the unresolved turbulence stresses r. or the scalar fluxes q t , will have to be 

provided by the subscale models discussed in the next section. It is not difficult to 
recognize that this second component will actually become the total Reynolds stresses 
and fluxes when the width of the temporal filter becomes sufficiently large, as the 
contributions from both T.. and T. will vanish. 
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1.4 Mathematical Models for Unresolved Turbulent Stresses and Fluxes 


In order to obtain the solution for the large scale turbulence using PRNS equations (4) - 
(6), we must close this set of equations by modeling the unclosed terms defined in 

Equation (7): z ij = 'p(u i u j -u i u : .), q i =p{u i e-u i e). They represent the effects of 
unresolved turbulence in the form of stresses and fluxes. 

There are many ways to model these unclosed terms. The more sophisticated method is to 
directly solve the transport equations of the unresolved turbulent stresses and fluxes, 
which can be derived from the Navier-Stokes equations (see Reference [8]). This method 
will require modeling of even higher order unclosed terms appearing in the quite complex 
transport equation system of r (/ and q l . A less complicated way is to start from a general 

constitutive relationship between the unresolved turbulent stresses r ;/ . and the strain rate 
of resolved large scale turbulence s tj , hr . This general constitutive relationship is then 

simplified according to the flow complexity by truncating the higher order nonlinear 
terms of s tj , hf . For example, the simplest form is just a linear relationship; and this is 

the widely used eddy viscosity model. Even at this level, there exists many approaches to 
provide the eddy viscosity. The simplest one is the Smagorinsky model [21] and its 
variations, which explicitly uses the grid size A as the length scale. A more sophisticated 
one is the one-equation model, such as the one proposed by Menon [19], which solves the 

transport equation of the unresolved turbulent kinetic energy k , and using yfk as the 
velocity scale, but still using the grid size as the length scale. In the PRNS approach, a 
two-equation model is adopted, which solves the transport equations of the turbulent 
kinetic energy k and its dissipation rate s , and uses them to form the length and velocity 
scales for the unresolved turbulence, thereby eliminating the use of grid size as the length 
scale for the unresolved turbulence. 

Here, we will present a constitutive relationship for the unresolved turbulent stresses, 
which is derived and simplified from the general constitutive relationship [22] by using 
the realizability and rapid distortion theory limit. This model contains linear, quadratic 
and cubic terms. 


1.4.1 Modeling of Unresolved Turbulent Stresses r 


The model proposed for PRNS is the following: 
T ij =- 2 fl C Mp^( S ij- 5 i S kk l 2 ) + \ S ijP k 

~ ^ 3 ./? P — ( k ik®kj ~ ®ik k kj ) 


+ 2 4/s P K [®ik~ s lj - ~ s l% + P>ik~ s kAj - <~ s iAAj + 11 s (s 9 - $ij~hk / 3 )] • 


(19) 
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Where, s y =(u IJ +u jJ )/ 2, a tJ =(u iJ -u J .)/2 , H s =(s kk s mm -s kl s lk )/2 . The model 
coefficients C , A 3 and A 5 are constrained by the realizability condition and rapid 
distortion theory limit. They are not arbitrary but formulated as (see Reference [22]): 


C ,=' 


1 


5 ^3 


1.0-A^I 

( k A 

— S 

J M 1 

U J 


4.0 + A S -U 
s 


0.5 + 1.5^ny 
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k 2 

1 - 6 C u p K - 


_ k 4 7S*S* + Q*Q* ’ 
P A 4 


in which, 


A'=j6cos(p, ^ = -^arccos(V6ff*j, 


r-MA 


(S) 


3 5 


+(«•)’- , 5-=^, £1-=^, 




( 20 ) 


( 21 ) 

( 22 ) 


The coefficients f x , / 3 and / 5 are functions of A y ./r , i.e. the ratio of the temporal filter 

width to the global time scale of the turbulent flow of interest. These functions must have 
the following property: 


A 



0 

1 


if 

if 



(23) 


This is because the unresolved turbulent stresses r y must vanish when the filter width A r 
vanishes, and r. must approach the Reynolds stresses R as A r increases towards T . In 
PRNS, A t /T< 1 , so we may make the following expansion: 


f, 


T j 


= c'o+ +c' 


aA 2 

T y 


+ • 


i = 1,3,5. 


(24) 


Where all C‘ 0 must be zero, because /’ must be zero as A r goes to zero. If we retain 
only the leading term as a first order approximation, then all will have the same form: 

f i (A T /T)*A T /T (25) 

We refer the ratio A t /T as the resolution control parameter (RCP). It controls the 
content of the resolved turbulent scales in the PRNS simulation. We will further discuss 
this parameter in Section 1.4.5. It should be pointed out that f. may not be all that simple 
as shown in Equation (25), more complex functional forms are possible, such as 
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(26) 


ft =[l-exp(-4'- 



A 


It satisfies the requirement of Equation (23), where A[, A l 2 are constants and need to be 
determined. 


1.4.2 Physical Effects of Unresolved Turbulence 

It is important to identify the various mechanisms occurring in the interactions between 
the resolved and unresolved turbulent scales, and how they are actually modeled in the 
numerical simulations. In the momentum Equation (4), these interactions are accounted 
for by the term z\. . , which is unclosed and must be modeled. In traditional LES, this 

term is modeled via the isotropic eddy viscosity. Therefore, the effect of unresolved, 
small scales appears to be only a modification to the viscosity of the fluid. However, the 
real, physical situation may be much more complex than just interaction via the eddy 
viscosity. In fact, the general constitutive relationship of r ;/ contains more than just the 

contribution from an eddy viscosity. As an example, Equation (19), a model for r. , 

contains two parts: linear term and nonlinear (quadratic and cubic) terms. Each part plays 
different role in the momentum equation. The linear part acts like additional viscosity 
(called subscale eddy viscosity); and the nonlinear terms act like additional source terms 
to promote the resolved large scale turbulence. This can be clearly identified by plugging 
T tj model into Equation (4), 


( 2 'l 

/ 

P + ~pk 

+ 

V 3 J 

J V 


2 

y 


1 ),j 


where, 


fr -k 2 

Pt =fx C »P — 


(28) 


Sf = \ A 3 f 3 p—(s^co u - co,,s, 


'ik^kj ^ik°kj 


.) 


(29) 


I —kr -i 

2 4/ 5 P — !_<%• S i - S ?k% + aJbn&mi ~ “>klV m Aj + 11 s / 3 ) J 


Apparently, the linear part of the model adds an additional subscale eddy viscosity ju T 
(which is isotropic) to the physical viscosity of the fluid /7 ; and the nonlinear part 
provides a complex source term Sf , which accounts for the effects of anisotropy and 
rotation. 
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We have noticed that, although different subscale eddy viscosity models have been used 
in different LES approaches, they all have neglected the source term Sf . Our preliminary 

study shows that the source term could become critically important for some flow 
simulations, especially at relatively low Reynolds numbers, and this will be briefly 
presented in Section 2.3. 


1.4.3 Subscale k- £ Transport Equations 

To complete the proposed model for t .. , we need k and s , the unresolved turbulent 

kinetic energy and its dissipation rate. Their exact transport equations can be derived 
from the Navier-Stokes equations, and contain some higher order unclosed terms due to 
the temporal filtering operation. Here, we briefly describe the procedure of the derivation. 
The first step is to establish the transport equation for r y , followed by tracing to establish 

the equation for z u (which is 2 pk), and this leads to the transport equation for k . 


The exact transport equation for the unresolved turbulent stresses zy. ( pu j n j - pu i u j ) is 


ij,< 


+ (^ T ij) k = D iJ + ®iJ +P ii-P £ ij’ 


(30) 


where Z> , ® ;/ , R and pe tj are the diffusion term, the pressure-strain term, the 

production term, and the dissipation term, respectively. The following expressions 
indicate that all terms on the right hand side of the equation, except for the production 
term R , are unclosed and must be modeled. 


A i = ~{p u , u j u k ~ P u , u j Rt ) + 1 2 P UjS ik - ^ 8 ik p u jSmm - 


~ 
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~( _ 2 _ ^ 

u j 2 Ps ik --S ik p Si 
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+ (r ik Uj + Tpu, ) k -(puj 8 ik + pu, 8 jk - p Uj 8 ik - p il, 8 Jk ) * 


f . 2 


( 2 \\) 

2 pu i s jk --8 jk pu i s mm - 

“i| 

2ps jk - — 8 jk p s mm 
v 3 J JJ 


® lJ =2ps ij -2pk ij 

k^j,k 

2p( 


P — —TrU. J —T, U J 

ij ik j,k jk i,k 


P £ ij 


s ik u j, k +s jk u i,k)--P s mm s ij 


2 p( 




Now, if the diffusion term Z> is modeled by a gradient-type diffusion of the eddy 
viscosity, then the trace of Equation (30) becomes 
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d d _ 3 /_ x 3 

—pk + —pu i k = —{p + p T )—k -rsy-ps 

ot ox,. Sr.. Sr.. 


in which, dx has also been neglected by ignoring the compressibility. 


The dissipation rate is defined by p s u / 2 , i.e., 


( _ 2 _ ) f _ 2 _ 

pe= 2 ps ij s ij --ps mm s ii - 2 ps ij s ij --ps mm s ii 
v j J v j 


A model transport equation for the dissipation rate s can be constructed by the analogy 
to Equation (31) as 


3 _ d __ S /_ „ s ps 

- ps + -pu,S = - -C mV ,j-C.,— 


where C £l and C e2 are the model coefficients. We have adopted the commonly used 
values of C sl = 1 .45 and C e2 = 1 .92 in the present work; keeping in mind that they can 
be further constructed as functions of local turbulence quantities [23]. 


1.4.4 Modeling of Unresolved Turbulent Fluxes q t 

A common practice in modeling the unresolved turbulent fluxes q t = p (up - up) is to 
employ the following isotropic model: 

q, = ~k t e t (34) 

k t is the eddy diffusivity for the heat, which is often modeled as k t = p T t Pr r , where 
Pr 7 (about 0.9) is the turbulent Prandtl number. However, based on the analysis of 

constitutive relationship (see References [24, 25]), the simplest form that considers the 
effects of strain and rotation should be 


q, = ~ K T e , ~ -(fAy + c 2 cOy )e j 


Where c x and c 2 are some undetermined coefficients. This more general model will add 
modifications to both the diffusion and source terms in Equation (5): 

r 2 _ 

(Pe) t + (P a i*) i = {(* + K T)h) - + P S * + 1 P s ij s i j --P s kk s n +Q + S t (36) 

V * J 
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where the extra source term from the unresolved turbulent fluxes is 


S e T = 


k( _ 

—\ c \ s ij 


T ^ y'V’ij C 2 f / 


( 37 ) 


So far, this term has not been considered in any of the CFD simulations that we are aware 
of, including the present work. 


1.4.5 Resolution Control Parameter RCP 


In the PRNS approach, the resolution control parameter RCP, i.e., the ratio A t /T , is 

used to regulate the scale content of the resolved field through its role in the models of 
the unresolved turbulent stresses and fluxes; it implicitly determines the low-end of the 
resolved, large scales. When RCP — » 1.0, the low-end of the resolved scales becomes the 
scale of turbulent mean flow field, therefore, the unresolved field covers all scales of the 
turbulence, and the PRNS simulation essentially becomes an unsteady RANS simulation. 
As the value of RCP decreases, the low-end of the resolved scales is expected to move 
towards small scales; and the unresolved field will only cover the rest of smaller scales. 

In planning a very large eddy simulation, we are often concerned with the following 
questions: how to choose the value of RCP from the outset followed by estimating the 
corresponding numerical grid size required for this simulation? Alternatively, for a given 
numerical grid size, how to choose RCP to yield a resolved scale level that the numerical 
grid can support? 


The following length scale definitions and relationships discussed in References [8], [26] 
and [27] may help to answer these questions. Let us denote k RANS , s R4NS and k, s the total 

turbulent kinetic energy, the total dissipation rate and the subscale turbulent kinetic 
energy, the subscale dissipation rate, respectively. A simple dimensional analysis 
suggests that an estimate for the smallest turbulent length scale that can be practically 

/ o / x 1/4 

resolved by PRNS is tj prns = \v T / s) , where v T is the subscale viscosity and s the 

subscale dissipation rate. In other words, rj PRNS represents the low-end of the resolved 
scales in PRNS. Similarly, the turbulent integral length scale is estimated to be of the 
order of £ MNS = 1/4 (this is consistent with a common estimation given by 

k 2 RANs/ £ RANs X where v RANS denotes the total turbulent eddy viscosity. Obviously, the 
mesh size A should be smaller than or equal to r/ PRNS , and the following relationship can 
be derived using the definition of RCP (see [8, 26, 27]): 


RCP> 


f A ] 

4/3 

r k ^ 

0 

V RANS J 


h 

\ * RANS J 


( 38 ) 


NAS A/TM— 2008-2 15418 


16 



When RCP is given from the outset, Equation (38) can be used to estimate the largest 
possible value of the grid spacing A relative to the turbulent integral length scale P RANS 

for a given ratio of k/k RANS (e.g., 10%). On the other hand, if the mesh size is given from 
the outset, Equation (38) can be used to estimate the value of RCP that the given grid size 
can support. For example, when a given grid resolution A! P RANS is about 1.5% and 

k/kRA NS about 11%, Equation (38) indicates that this grid resolution can support PRNS 

simulations having RCP from 1.0 to a value as small as 0.3. For the same given grid, we 
would obtain RANS type results if RCP were about 1.0, and LES or VLES type results 
if RCP were about 0.3 - 0.4. This has been demonstrated in the simulation of the flow 
field in a single injector flame tube (see Section 3.1). The vorticity contour plots clearly 
illustrate that the simulated flow field transits from RANS type toward LES type when 
RCP varies from 1.0 to 0.38. 

Furthermore, let L be the length of the computational domain, A rans be the mesh size 
used in the RANS simulation, then the number of grid points needed for a PRNS 
simulation, N prns ~ (f/A) 3 , can be estimated as 

AW > AW ' (A„« /(„ ) ! ■ Rep* (t/t„ Y ,n , (39) 

where « (Z/A^^) 3 represents the number of grid points required by RANS. 

Equation (39) provides an estimation of the number of grid points needed for a PRNS 
simulation when the ratio A RANS j P RANS (e.g., 2%) and k/k RANS (e.g., 10%) are given. 


2. Evaluation of PRNS using NCC 

Ultimately, we want to use the NCC for performing the very large eddy simulation of 
combustor flows. At the present stage, the focus is the fundamental assessment of the 
implemented subscale model and the establishment of a sufficiently “clean” numerical 
setting for the PRNS calculations. 

We have conducted systematic studies to examine the effects of the following numerical 
artifacts on the calculation of turbulent fluctuations: the truncation errors associated with 
the temporal-spatial discretization, the convergence errors associated with the “pseudo- 
time” iteration, the generation of disturbances by forcing the redistribution of the 
residuals, the round off errors associated with the computation of the pressure gradient in 
low Mach number compressible flow, and the added artificial dissipation required for 
maintaining the numerical stability. The results of these studies indicate that the 
calculated turbulent fluctuations are quite sensitive to the setting of these numerical 
artifacts, and we have established a sufficiently “clean” numerical setting for the PRNS 
calculations (see the following section). The simulations presented in this paper are all 
conducted using this numerical setting as the baseline. 
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At this stage of the evaluation, we have selected the fully developed turbulent pipe flows 
as the focus, because they are well known fundamental flows with experimental data [20], 
and they are amenable to the use of periodic boundary conditions, thereby, allowing us to 
concentrate on the subscale model itself and free from the complicating issues associated 
with the inflow and outflow boundary conditions. We will at first examine the PRNS 
results obtained from using only the eddy viscosity type of model at different Reynolds 
numbers, followed by exploring the behavior of a more general subscale model that 
accounts for both the effect of eddy viscosity and the effects of anisotropy and rotation. 

In this report, we will present the results obtained using various values of RCP for a 
moderately large Reynolds number of 150,000 and a relatively low Reynolds number of 
15,000 (based on the centerline axial velocity and the pipe diameter). Our criteria for a 
successful simulation are: the simulation must produce a broadband scales of turbulent 
fluctuations, in particular, the power spectrum must cover at least two or three orders 
energy variations from large to small scales; and the simulation must produce reasonable 
or better mean profdes than those provided by the RANS simulation, when compared 
with available experimental data. 

In the following, Sections 2.1 and 2.2, discuss the PRNS results obtained by using a 
subscale model accounting for only the effect of eddy viscosity, i.e., the unresolved 
turbulent stress x tj is represented only by the linear part of Equation (19). Subsequently, 

the results of a preliminary study are described in Section 2.3, which demonstrate that the 
nonlinear part of r ;/ can have a significant effect on the numerical solution, especially for 

flows at lower Reynolds numbers. 


2.1 Turbulent Pipe Flow at Reynolds Number 150,000 

A fully developed turbulent pipe flow at Reynolds number 150,000 has been selected to 
represent turbulent pipe flows at moderate or large Reynolds numbers, in which the 
turbulent fluctuations are quite energetic and will not be easily damped by the inherent 
numerical dissipation of a CFD code. Figure 1 shows the computational domain of the 
pipe with a diameter of 0.12936 m, which has an aspect ratio (length to diameter) of five 
and 900,705 hexahedra elements. Such a grid resolution is typical for a fine grid RANS 
simulation. We will use this same grid for all PRNS simulations. 
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Figurel: Computational domain for the pipe flow simulations. 

Periodic boundary conditions are imposed at the inlet and the outlet of the computational 
domain. A generalized wall function [17] is applied for the solid wall. The initial 
conditions are created by a specified axial velocity profile randomly disturbed in all three 
directions. All simulations are performed using the NCC code, which uses the dual time 
stepping method. The baseline numerical setting is the following: 

• Initial gauge (fluctuating) pressure = 0.0. 

• No residual smoothing is applied. 

• Coefficients of the second- and fourth-difference dissipation terms (2 nd and 4 th ) 
are 0.0 and 0.001, respectively. 

• For the convergence of the pseudo-time (inner) iteration, three orders of 
magnitude reduction in the residuals is typically required, i.e., the convergence 
order = 3. Otherwise, it is stipulated that the iteration must reach the maximum 
number of 120. The CFL number for the pseudo time is 1.0. 

• The physical time step dt = 2.5E-05 second. 

Following the analysis described in Section 1.4.5, in order for the present grid resolution 
to be able to capture those turbulent fluctuations that are responsible for more than 90% 
of the total turbulent kinetic energy, the value of RCP should be between 0. 3-0.4. If RCP 
increases, the content of the unresolved scales will expand into larger eddies, and the 
unsteadiness of the calculated, resolved part will diminish. Keep increasing RCP, the 
PRNS solution will eventually become the RANS solution. This trend can be observed 
from the numerical results with RCP = 0.2, 0.3, 0.34, 0.4, 0.5, 0.8 and 1.0, which are 
presented in sections A.1-A.7 of the Appendix A. Furthermore, the results obtained by 
the traditional RANS approach are also shown in section A. 8. Note that the numerical 
setting for RANS is somewhat different from the numerical setting of PRNS, but the 
same grid is used. 
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Here, we present the numerical simulation with RCP = 0.3 as a representative of PRNS 
simulation of very large scale turbulence. We consider simulations with RCP > 0.4 as 
the very, very large eddy simulations. 


Figure 2 shows the time histories of velocity components recorded at the probe 1, which 
is located at the centerline of the pipe. Here, u and v are the tangential and the radial 
velocity components, respectively; and w is the axial velocity component. Time histories 
at other locations, such as the probe 6 (the point between the centerline and the wall) and 
the probe 14 (near the wall), can be found in Section A.2. These time histories indicate 
that the turbulent fluctuations are fully developed after 10,000 time steps. The fully 
developed, various turbulent scales can also be seen from the snap shots of various 
turbulent quantities, such as u, v, w, pg, k, and ju + /u T . For example, Figure 3 is the 

contour plots of w and u at time step 36,400, and Figure 4 is the contour plots of gauge 
pressure pg and subscale turbulent kinetic energy k. These snap shots do reveal the 
presence of wide range turbulent scales and the statistical homogeneity in the axial 
direction. As a fully developed pipe flow, the tangential and the radial velocity 
components must be statistically zero, and the gauge pressure (which does not include the 
mean pressure gradient) must be statistically constant. These conditions are in deed met 
by the numerical simulations, as shown by both the time histories and the contour plots. 
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Figure 2: Time histories of velocity components at the centerline (Probe 1). 
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Figure 3: Snapshot contours of w and u at the time step 36400. 
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Figure 4: Snapshot contours of gauge pressure and subscale turbulent kinetic energy. 


Examining the time power spectrum of turbulent fluctuations and the two-point 
correlations is another way to assess the quality of the numerical simulation. Figure 5 and 
Figure 6 present the power spectrum density (PSD) of the w component and its two-point 
(time) correlation at three locations (Probes 1, 6 and 14). The broadband feature of the 
PSD (i.e. more than two orders of energy variation from small scale to large scale) and 
the typical two-point correlation shapes (i.e. the correlation rapidly decreases as the time 
lag increases) indicate that the PRNS simulation with RCP = 0.3 does mimic the 
statistical features of a fully developed turbulence. 
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Figure 5: Power spectrum density of w component at Probes 1, 6 and 14. 




Figure 6: Two point (time) correlations of <ww> at Probes 1, 6 and 14. 


There exists an experimental data of axial mean velocity profile at Reynolds number of 
145,700, which is very close to 150,000. (Reference [20]). Figure 7 compares the PRNS 
results with the experimental data and the RANS results. As it can be seen, the result of 
PRNS compares reasonably well with the experimental data while the steady RANS 
result exhibits significant under-prediction. 
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Pipe flow, Re = 150,000 


0 PRNS with RCP = 0.3 



Figure 7: Comparison of W mean profiles between PRNS, RANS and experimental data. 


To examine the behavior of PRNS at different values of RCP, we have also performed 
simulations with RCP = 0.2, 0.34, 0.4, 0.5, 0.8 and 1.0 using the same numerical setting. 
The results of these individual cases can be found in Appendix A. In the following, we 
summarize the major observations distilled from this parametric investigation. 

In general, we find that, for a given grid resolution, there exists a minimum value of RCP 
(say 0.3) below which (say 0.2) the intended scale content of the operating subscale 
model cannot be properly supported by the resolution level of the given grid. In effect, 
this amounts to an attempt of using a coarse grid to execute a direct numerical simulation 
(DNS). The results will not be physically meaningful, and this is demonstrated by the 
results of the simulation with RCP = 0.2 (Section A.l). Although we still see a strong 
development of fluctuations, the operating subscale model does not play a role in their 
development, because the subscale turbulent kinetic energy appears to be erroneously 
small, i.e. k is of order l.E-1 1, and the subscale eddy viscosity is nearly zero. 

On the other hand, when RCP increases (e.g., RCP = 0.3, 0.34, 0.4... 1.0), the resolved 
part will change from fluctuations consisting of very large scale turbulence toward 
Reynolds-averaged mean value. This can be seen from Figure 8, the time history of the 
centerline axial velocity component (w), and Figure 9, the snapshot of the w contours at 
the end of each simulation. Figure 10 shows the instantaneous axial velocity w profdes 
from PRNS simulations with RCP = 0.3, 0.5, 0.8 and 1.0, also included is the 
experimental data for the mean profile. These figures indicate that PRNS with RCP = 0.3 
can mimic the turbulent motion associated with a wide range of large scale turbulence, 
and can predict the mean profile reasonably well. 
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Figure 9: Contours of w component at center plane with RCP = 0.3, 0.5 and 1.0. 
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Figure 10: Instantaneous w profiles with RCP = 0.3, 0.5, 0.8 and 1.0. 
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The results shown in Figure 10 also indicate that the accuracy of the current subscale 
model becomes progressively worse as RCP increases, and, when RCP=1, the time- 
averaged profile coincides with the mean profile from the steady RANS calculation using 
a standard k-s model (see Section A.8). This is expected, because, when RCP = 1, the 
current PRNS is equivalent to an unsteady RANS simulation using a standard k-s 
model. It is desirable to further improve the current PRNS subscale model so that it will 
lead to better RANS results than those from the standard k-s mode. For example, a so 
called “Low-Reynolds number k-s model” ([28]) may be considered as a candidate for 
this purpose. 


2.2 Turbulent Pipe Flow at Reynolds number 15,000 

A pipe flow at Reynolds number 15,000 has been selected to represent turbulent pipe 
flows at lower Reynolds numbers. In these cases, physically, the turbulent fluctuations 
are relatively weak, and computationally, they are more susceptible to be damped by the 
numerical dissipation associated with the truncation errors and the added artificial 
dissipation. We have performed the PRNS simulations using RCP = 0.0, 0.2, 0.3, 0.34, 
0.4, 0.5 and 1.0. The grid and the numerical setting are the same as those for Re= 150,000, 
except that the time step is now dt = 4.E-04 (the flow is now slower) and a much smaller 
added fourth-difference dissipation (4 th = 0.00001) is used to further reduce the numerical 
damping. The detailed, individual results can be found in Appendix B (B.1-B.7). 
Generally speaking, in all of these simulations, the development of the initially induced 
fluctuations can not be sustained and the calculated flows do not evolve towards a fully 
developed turbulent flow. 

In the case of RCP = 0.0, the simulation is performed without the use of subscale eddy 
viscosity. The results in Appendix B.l indicate that the turbulent fluctuations have been 
damped out and the flow is evolving toward a laminar pipe flow. 

In the case of RCP = 0.2, we see that the initially induced fluctuations are damped out 
sooner than in the case of RCP =0.0. This is apparently due to the effect of the invoked 
subscale eddy viscosity. As shown in Appendix B.2, all velocity fluctuations are damped 
out, and the subscale turbulent kinetic energy is practically zero. Therefore, the 
simulation is also evolving toward a laminar pipe flow. 

When RCP = 0.3 or larger, the simulations appear to be quite different. The subscale 
turbulent kinetic energy is no longer vanishingly small, and the operating subscale eddy 
viscosity is more appreciable. 

In the case of RCP = 0.3 (Appendix B.3), we see some sustained fluctuations. But clearly, 
it is not evolving toward a fully developed turbulent pipe flow. It seems that, at this 
particular value of RCP, the operating subscale eddy viscosity is small enough to still 
allow the presence of some large scale turbulence in the resolved field, but large enough 
to hinder their further development. 


NAS A/TM— 2008-2 15418 


25 



In the cases of RCP = 0.34, 0.4 and 0.5 (Appendix B.4, B.5, and B.6), we see that all 
fluctuations are again damped out, but the unresolved turbulent kinetic energy remains 
appreciably large. This suggests that the calculated solution is the “mean” of some under- 
developed turbulent flow, in which the dissipation effect of the operating subscale 
viscosity becomes unrealistically dominant. 

In the case of RCP =1.0 (Appendix B.7), the PRNS simulation amounts to an unsteady 
RNAS simulation. The calculated solution is the Reynolds-averaged mean value. The 
accuracy of the PRNS results is the same as that of the steady RANS results using a 
standard k-e model. 

It should be noted here that this Re= 15,000 pipe flow had been previously simulated by 
PRNS using the NCC with RCP = 0.38. Although the same grid was used, the numerical 
setting was quite different, namely, the previous setting was not as computationally 
“clean” as the present one. It had allowed the round off error associated with the 
computation of the pressure gradient to become appreciable, it had adopted a lower order 
convergence criterion for the inner “pseudo-time” iteration, and it had used residual 
redistribution during the inner iteration. The net effect was that various, apparently 
random sources were injected into the simulation to promote the appearance of 
fluctuations in the calculated results. Somehow, the time-averaged results compared 
reasonably well with the experimental data (see Reference [8]), we now consider the 
good comparison as fortuitous. Nevertheless, this experience does suggest that, in 
addition to the eddy viscosity type of effect, physics-based sources type of effect need to 
be included in the subscale model to more faithfully mimic the interactions between the 
resolved and un-resolved turbulence scales, in particular, when the physical fluctuations 
involved are relatively weak, or when the local Reynolds number is relatively low. We 
have been exploring this idea, and some of the preliminary results are presented in the 
following section. 


2.3 Preliminary Study of A Nonlinear Subscale Model 

In Section 1 .4, we have provided a nonlinear formulation of the subscale model for the 
unresolved turbulent stresses r. , Equation (19), as well as the resulting source terms in 

the momentum equation and the energy equation, i.e., Equation (29) and Equation (37), 
respectively. Equation (19) may be rewritten as 

hj - 1 #ijP k = - 2 Pt (s v - SyS kk /3) 

— ^3 A P ~~2\^ik^kj — ^ik^kj ) (40) 

+ 2v4 5 / 5 p —\jn ik s kj — s ik (b kj + cb ik s km S) mj — o) kl ^ lrr fi> mk S^ + II s (Sy — SySjx /3)J . 


NAS A/TM— 2008-2 15418 


26 



where // 7 is the subscale eddy viscosity defined by Equation (28). And, at this stage of 
the development, all the coefficients f x , / 3 and f 5 are set equal to the resolution control 
parameter RCP. 

The PRNS results presented in Sections 2.1 and 2.2 are obtained with a subscale model 
involving only /u T , the nonlinear part that reflects the effects of anisotropy and rotation 

has been ignored. As mentioned before, this nonlinear part leads to additional source term 
in the momentum equation for the resolved turbulence, and this source term may play an 
important role when the unresolved scales become larger or when the flows are at lower 
Reynolds numbers. 

To explore this idea, we have included the nonlinear terms in the subscale model and 
repeated the PRNS simulation of the pipe flow at Reynolds number 15,000. The results 
are very encouraging. The detailed results for the case of RCP = 0.3 can be found in the 
Appendix C. 1 . They indicate that the fluctuations of velocity components are now being 
sustained and evolving towards a fully developed state. This is clearly due to the new 
source terms from the nonlinear part of r j; . 

The effect of nonlinear subscale model on the PRNS simulations for flows at relatively 
large Reynolds numbers has been examined by using the pipe flow at Reynolds number 
150,000 with RCP = 0.3. For this Reynolds number, we employed the existing flow field 
previously established with the linear subscale model as the starting condition, and the 
simulation was then performed with the nonlinear subscale model. We have not found 
much dramatic changes between the two simulations, except that the subscale turbulent 
kinetic energy as well as the subscale eddy viscosity appears to be stronger in some 
regions. The detailed results can be founded in Appendix C.2. 

A more extensive study of the effects of the nonlinear subscale model is under way, and 
the results will be reported elsewhere. 


3 Preliminary Applications 

To illustrate the intended applications for the PRNS approach, we include here some 
sample results previously published in References [8] and [27]. 


3.1 Non-Reacting Flow in A Single Injector Flame Tube 

This injector has been used in an industrial combustor. The flow field inside a single 
injector flame tube has many important features occurring inside a real combustor, such 
as massive separation, strong swirling and recirculation. This flow has been studied both 
experimentally and numerically [29], 
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Figure 1 1 : Computational domain for RANS and PRNS simulations. 

Figure 1 1 shows the computational domain. The combustion chamber is a rectangular 
box. A highly swirling jet is injected from a circular inlet. The Reynolds number based 
on the inlet axial velocity and diameter is about 3,200,000. A grid of 495,000 elements 
was used in the simulation. At the inlet, only the mean profdes of velocities, density and 
temperature are specified. At the outlet, simple extrapolation is used to minimize the 
possible unphysical reflections. Figures 12 and 13 are the instantaneous contours of the 
axial velocity and the vorticity magnitude, respectively. They are obtained from three 
different values of Rep : 1.0, 0.46, and 0.38. Figures 14 and 15 illustrate the 
instantaneous features of the kinetic energy of the unresolved scales and the effective 
eddy viscosity (p + p T ). These results demonstrate that the characteristics of simulation 

will change from RANS towards LES as the RCP reduces from 1.0. Figure 16 gives the 
comparison between several instantaneous centerline axial velocity distributions with the 
mean value provided by LES data. Apparently, the simulation with RCP = 0.38 and a 
relatively coarse grid captured many features of LES. Figure 17 shows the time history of 
the axial velocity fluctuation and its power spectrum at a centerline location inside the 
recirculation region. The broadband feature of the power spectrum clearly indicates that a 
wide range of large-scale turbulence has been captured by PRNS. 
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Figure 12: Contours of axial velocity Figure 13: Contours of vorticity magnitude. 
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Figure 14: Unresolved kinetic energy. 


Figure 15: Subscale eddy viscosity. 


NAS A/TM— 2008-2 15418 


29 



150 - 


PRNS with 2-Eq. Rcp=0.38 



Figure 16: Axial velocity profiles. Figure 17: Velocities and power spectrum. 


3.2 Reacting Flow in a Lean Direct Injection (LDI) Hydrogen Combustor 

The so called C4 configuration was one of the experimental LDI hydrogen combustors 
[30] studied at NASA GRC. It was designed to study the effects of fuel/air mixing on the 
combustion, especially on the NOx emission. In general, better mixing of fuel and air will 
improve the combustion as well as reduce the NOx emission. To assist the experimental 
studies, we had carried out numerical simulations using both RANS and PRNS 
approaches. 

The geometry for the C4 configuration is shown in Figure 18. It consists of seven 
converging-diverging nozzles, the fuel (hydrogen) is injected into the upper part of each 
nozzle perpendicularly through four fine fuel tubes, the air from the inlet flows to nozzles 
and meets with the high speed hydrogen jets, the mixture of air and hydrogen then flows 
through the converging-diverging nozzles into the combustor. The studies were focused 
on the effectiveness of fuel/air mixing and its effect on combustion. An unstructured 
tetrahedral grid was used for this complex geometry. The total number of cells is about 
3.5 millions. The same grid has been used for both RANS and PRNS simulations. 



Figure 18: The configuration of a LDI hydrogen combustor. 
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RANS Simulation of the C4 LDI Hydrogen Combustor 

We started with the steady RANS simulation of a non-reacting flow using the nonlinear 
k-s RANS model [15]. The aim was to examine the global features of the fuel/air mixing. 
Subsequently, we switched the combustion on to perform the steady RANS simulation of 
the reacting flow. The combustion model used in this simulation was a one based on the 
eddy dissipation concept [31]. The inlet and outlet boundary conditions were set 
according to the experimental studies. The mass flow rate and temperature for air at the 
inlet were 0.40012 kg/s and 700.212 K; the mass flow rate and temperature for hydrogen 
at its inlet were 0.00482 kg/s and 588.555 K. The outlet back pressure was set to 622,040 
Pa. 


Results of steady non-reacting RANS simulation 


The global flow feature is indicated in Figure 19, which is a contour plot of the axial 
velocity at a center plane. Figure 20 shows the vorticity magnitude at the same plane. 
Figure 21 is the distribution of the hydrogen mass fraction in the lateral cross section at 
the nozzle throat. Figure 22 is the distribution of the hydrogen mass fraction at a different 
lateral section located at the exit of the nozzles, i.e., the entrance of the combustor. These 
fuel mass fraction plots can be used to assess how well the mixing is achieved by the 
design of this particular C4 configuration. 
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Figure 19: Axial velocity (RANS). 


Figure 20: Vorticity Magnitude 
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Figure 21 : Fuel mass fraction at nozzle throat. Figure 22: Fuel mass fraction at nozzle 

exit. 


Results of steady reacting RANS simulation 


As mentioned before, a simple turbulent combustion model, namely, the Magnussen 
model [31] was used to obtain a global picture of the flame structures. Figures 23 and 24 
are the contour plots of axial velocity and vorticity magnitude at the center plane. We 
observed that the axial velocity and the vorticity are much larger and stronger than that 
revealed in the non-reacting simulation. Figure 25 is the distribution of temperature at the 
center plane, and Figure 26 is the temperature contour in the neighborhood of a nozzle. 
These temperature distributions indicate that the combustion occurs quite close to the exit 
of the nozzles. 
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Figure 23: Contour of axial velocity. 


Figure 24: Contour of vorticity magnitude. 
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Figure 25: Temperature at center plane. 


Figure 26: Temperature near the injection 


PRNS Simulation of the C4 LDI Hydrogen Combustor 

A very large eddy simulation using PRNS was carried out using the steady RANS results 
of the reacting flow as the initial condition. The resolution control parameter was RCP = 
0.38. 

Figures 27 and 28 are the snapshots of the contours of the axial velocity and the vorticity 
magnitude, which indicate the presence of unsteady, large scale flow structures that are 
absent from the corresponding RANS simulation (see Figures 23 and 24). Figure 29 
shows the instantaneous temperature contours at the center plane. Figure 30 is the 
temperature distribution in a nozzle. These figures also show a very different flame 
structures from that provided by the RANS simulation (see Figures 28 and 29). We also 
noticed that the maximum flame temperature, 2543 K, as predicted by PRNS near the exit 
of the nozzles, is lower over 200 degrees than that predicted by RANS. 
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Figure 27: Contour of axial velocity. 


Figure 28: Contour of vorticity magnitude. 
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Figure 29: Temperature at the center plane Figure 30: Temperature near the injection. 


The RANS and PRNS simulations yield different features. RANS can quickly provide a 
global picture of the combusting flow, and can be used as the initial condition for the 
PRNS simulation. The explicit presence of the unsteady, large flow structures in the 
calculation will then improve the fidelity of the combustion simulation. By leveraging the 
strengths of these two complimentary approaches, higher fidelity analysis of turbulent 
combustion in realistic combustor and operating conditions can then be routinely 
conducted at reasonable computational costs. 


III. Findings and Recommendations 

The partially resolved numerical simulation (PRNS) approach for the very large eddy 
simulation (VLES) has been implemented into the National Combustion Code (NCC). 
Under the present effort, the fundamentals of the PRNS have been assessed by 
systematically examining the computed results of the fully developed pipe flows at 
different Reynolds numbers. The benefits of the PRNS for computing the practical 
combustor flows have been demonstrated by two representative cases. The major findings 
and recommendations emerging from the current investigation are: 

1. Major Findings 


• PRNS is effective in bridging the fidelity and the computing resource gaps between 
the traditional RANS and LES approaches. This feature should be particularly 
attractive to the researchers and designers working in an engineering environment. So 
far, our experience suggests that, with care, typical RANS type grid resolution can 
support PRNS to successfully perform the very large eddy simulation. 
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• PRNS has been verified to be a unifying approach for simulating the turbulent flows. 
It enables the computations evolving from the RANS toward YLES by controlling the 
value of the resolution parameter RCP. 

• Successful VLES simulations have been obtained with RCP between 0.3 and 0.4 for a 
very wide range of the flow Reynolds numbers. 

• Numerical setting in the NCC code, in fact, in any CFD code, must be “clean” enough 
such that the effects of the operating physical model are not appreciably polluted by 
the effects of the numerical artifacts. We have found a baseline setting for NCC to 
perform PRNS computations. 

• In lower Reynolds number flows, the physical, turbulent fluctuations are relatively 
weak, and they are more susceptible to the damping effects of both the numerical 
errors and the eddy viscosity of the subscale model. Additional physics-based sources 
are helpful in computationally sustaining the development of these weak fluctuations 
in the PRNS calculations. 

• A nonlinear formulation of the subscale model for unresolved turbulent stresses has 
been proposed. Its linear part accounts for the dissipation effect of the unresolved 
turbulence via the eddy viscosity, and its nonlinear part reflects the effects of 
anisotropy and rotation. The results of a preliminary evaluation indicate that the 
nonlinear part can play a significant role in sustaining the large scale turbulence, and 
its importance becomes more apparent in the case of lower Reynolds number flow. 

• Current RANS with standard k-s model in the NCC tends to under predict the 
mean axial velocity of the fully developed pipe flow. 

2. Recommendations 

• Comprehensive evaluation of the current nonlinear subscale model and its further 
development are needed to buttress the PRNS approach for the VLES of practical 
combustor flows. 

• A better RANS model is needed as the parent of the PRNS subscale model to 
improve the prediction capability of PRNS when the value of RCP approaches one. 

• Systematic study of the effects of grid quality and grid resolution on the PRNS 
simulation must be conducted. 

• Development of turbulent combustion models which are consistent with the scale 
resolution level of the PRNS approach is required for the predictive simulation of 
advanced combustion systems using the PRNS. 
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U, V 


Appendix A: Effects of the Resolution Control Parameter on the PRNS 
Solutions (Re= 150,000) 

A.l PRNS with RCP = 0.2 


History of velocity components at probe 1 and probe 14: 



Time Steps 


Time Steps 



Time Steps 


Time Steps 
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History of subscale turbulent kinetic energy at probes 1 and 14. 


IE-10 : 
SE-11 : 
SE-11 : 
7E-11 : 
6E-11 - 

5E-11 : 

* 

4E-11 1 
3E-11 | 
2E-11 : 
1IE-11 1 

0 : 

-IE-11 


Pipe Flow, Re= 150000, £ Probe 1 



10000 15000 

Time Steps 


IE-10 
0E-11 
0E-11 
7E-11 
5E-11 
5E-1 1 

'XL 

4E-11 
3&11 
2E-1 1 
IE-11 

0 

-IE-11 


Pipe Flow. Re= 150000, st Probe 14 
PRNS, RCP=0.2 

2nd=0.0, 4th=Q.001, CFL=1, dt=2.5e-5 
Co iwergen ce order= 3, or M ax iter= 1 20 




5000 10000 15000 20000 25000 

Time Steps 


Contours of w, u, v, pg, k, ju + p T at the time step 25,000: 


PRNS, Rep =0.2 

Re= 150000, 25000 time steps 

NCO-1 .1 A. no smoothing 

P-referen ce= 1 00290 

2nc=0.0, 4th=D.Q01 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



W 

2 49006+01 

2.331P6+Q1 
2 i923E+oi 
1.99358+01 
l 02408+01 
1 .05908*01 
1. 40738+01 

1 . 31006+01 

1 1400E+O1 
y 01Q2L>00 
6 1227E+QQ 

1 0 43516+00 
4.74706+00 
3.0001 e+oo 


PRNS, Rep =0.2 

Re= 150000, 25000 time steps 

NCO-1 .1 A. no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=D.Q01 p CFL=1.0, dt=2.5e-5 

Convergence 3 order or Max tter=120 



1 .70108+00 
1 47226+00 
1 10258*00 
Q 533 7 £-01 
5 43976-01 


2.34598-01 
-7 40476.02 
-3.W25&01 
a gaaarni 
-1.00318*00 
-1 312SE+Q0 


I 


-1 0219L*QO 
-1 .03136+00 
-2.24076+00 



PRNlS, Rep =0.2 

Re= 150000; 25000 time steps 

NCO-1 .1 A. no smoothing 

P-referen ce= 1 00290 

2nc=0 0, 4th=D.001 p CFL=1.0, dt=2.5e-5 

Convergence 3 order or Max tter=120 



1.54018+00 
1 .29046+00 
1 03056+00 
i 75U9L-D1 
5 .21156-01 
2.09416-01 
1 1070602 
•2 43078,01 
4 9701801 
-7.52556-01 
-1 .0073 E+QQ 


I 


-1 .29208+00 
-1 .0100E+OD 
.1.77158+00 


PRNS, Rep =0.2 

Re= 150000, 25000 time steps 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 00290 

2nc=Q.O, 4th=D.0G1 p CFL=1 .0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



I 


n 

7. £4348+00 
G. 45996+00 
E 37556+00 
4.293UL'UU 
3 2039E+0O 
2 12918+00 
1 04276+00 
4.09008 02 
1 1242r>00 
-2.20776+00 
-3 291 1 E+QQ 
-4 37406+00 
-5 40008+00 
-G £4 108+00 
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PRNS, Rep =0.2 

Re= 150000, 25000 time steps 

NCO-1 .1 .4, no smoothing 

P-referen ce= 1 00290 

2nc=0 .0 , 4th=Q.0G1 , C F L= ' 1 . 0, dt=2 . 5e-5 

Convergence 3 order or Max iter=120 



I 


2.0234&1 1 
2.03Q1&1 1 
2J370e.tl 
2.MS9&T 1 
2 0533601 
1. 0011&41 

i eases 1 1 
14705&11 
1 2043F 1 1 
i ua^uL-i i 
8 .0973 E- 1 2 
7.DM9&12 
S 1S20E-1 2 
£22B9&13 


PRNS, Rep =0.2 

Re=1 50000 , 25000 time steps 

NCO-1 .1 .4, no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=Q.0O1 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



I 


mu 

1.6970 SOS 
1.9970 EOS 
1 0970 £ 05 
1 U47UL-tf9 
1 007OE-O5 
1 WUL-0i 

i.satn&os 

1 .907O&Q5 
I fl07nr ns 
1.B970&D5 
1 897Q E-05 
1 UlrUL-lfi 
1 .6870645 
1.007Q&O5 



A.2 PRNS with RCP = 0.3 


History of velocity components at probe 1 and probe 14: 


Pipe Flew. Re= 150000, at Probe 1 
PRNS, RCP=0.3G 

2nrf=0O, 4th=0.001, CFL=1, dt=2.5e-5 
Convergence order= 3, or Maxiter=12G 



Pipe Flow. Re= 150000, at Probe 1 
PRNS, RCP=0.30 

2nd=0.Q, 4th= 0.001, CFL=1, dt=2.5e-5 
Convergence order= 3, or Maxiter=12Q 



20000 

Time Steps 


Pipe Flow, Re= 150000, d Probe 14 
PRNS, RCP=0.3G 

2nd=0.0, 4th=O.0O1, CFL=1. dt^2 5e-S 
Co nvergen ce order= 3. or M ax ifier= 1 20 



10000 20000 
Time Steps 


Pipe Flow, Re= 150OOO, d Probe 14 
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History of subscale turbulent kinetic energy at probes 1 and 14. 


0.16 r 


0.14 - 


Pipe Flow. Re= 150000, at Probe 1 
PRNS, RCP=0.3G 

2nd=Q.O, 4th= 0.001, CFL=1. dt=2.5e-5 
Convergence order= 3, or Max iter=120 


0.12 7 
0.1 r 
^ 0 .0 8 - 
0.06 7 
0 .04 +- 
0.02 - 


h 


3 

40000 



Pipe Flow. Re-150000, at Probe 14 
fi r PRIMS, RCP=0.30 

2nd=0.Q, 4th= 0 .00 1 , CFL= 1 . dt=2. 5e-5 
Co iwergen ce order =3, or M ax iter= 1 20 



0 10000 20000 30000 40000 


Time Steps 


Time Steps 


Contours of w, u, v, pg, k, // + ju T at the time step 36,400: 


PRNS, Rep =0.3 

Re=1 50000 p 36400 time steps 

MCO-t .1 .4. no smoothing 

P-referen ce= 1 0D29Q 

2nc=0.0, 4th=D.001 p CFL=1.0, dt=Z5e-5 

Convergence 3 order or Max iter=120 



i 

i 


2.33176*01 
2.17336+01 
2 0140ET + D1 
1 .05056+01 
1 .89916*01 
1.53976*01 
1.30136*01 
1.22296+01 
1 0G45r+Dl 
y 00126+UU 
7.47:26*00 
5.03336*00 
4 30346+0D 
2.72546+00 


PR MS, Rep =0.3 

R e= 1 50000 f 35400 time steps 

MCO-1 .1 .4. no smoothing 

P-referen ce= 1 0D29G 

2nc=0.0, 4th=D.001 p CFL“1.0. dt=Z5e-5 

Convergence 3 order or Max iter=120 



I 


i 


2.14056*00 
1 .72956+00 
1 3155E+00 
9 LW40L-(?1 
4 &2456P1 
IJ.LW4Ut.-5J2 
3 .3 1576.01 
7 4359601 
1 1 WISE* 00 
A .50706*00 
-1 .3709 E+QQ 
-2 .331 36*00 
•2.00396+00 
■3J1556+00 



PR MS, Rep =0.3 

Re=1 50000 p 35400 time steps 

MCC-1 .1 .4. no smoothing 

P-referen ce= 1 OD29Q 

2nc=0.0, 4th=0.001 p CFL=1.0, dt=2.5e-5 

Convergence 3 order or Max iter=120 



X 

2 10346*00 
1.0390 6+00 
i si07rHDn 
1 13306+00 


u roaa^.ai 


■ &4012L-O1 
2 2529 E-G1 

■ $ 70076-02 
w ' A 7047^01 

■7 4333 1-31 
-1 .0902 E+OQ 

U m -1 .30916+00 
-1.71196*00 
.2.0349 6+00 


PRMS, Rep =0.3 

Re=1 50000 p 35400 time steps 

MCC-1 .1 .4. no smoothing 

P-referen ce= 1 0D29G 

2nc=0.0, 4th=D.001 p CFL-1.0. dt=Z5e-5 

Convergence 3 order or Max tter=120 



I 


Pg 

9 30596*00 
7 29456+00 
9 20316+00 
5.27176+00 


4 I'srni njii 


■ 3.24026*00 

■ 2 23756*00 

■ 1 2291 E+QQ 
2 14546-01 
■7.2075601 
-1 .00026+00 

I ■ -2.01996*00 
H -3 8310E+00 

■ 4 04246+00 
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PR NS, Rep =0.3 

Re= 150000, 35400 time steps 

NCO-1 .1 .4. no smoothing 

P-referen ce= 1 00290 

2nc=0.0, 4th=0.001, CFL“1.0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



I 


lurtv 

i.doooe+oo 

e.23S5&Qi 

bjtwe-di 

7.71S4&01 
r gsae-E-m 

0 1323&01 
S43oe&oi 
4.0092&O1 
3 00TT&O1 
3 1402&D1 
2 3M0&O1 

1 0231S-31 
6.01S4&O2 
1 0000 E 02 


PR NS, Rep =0.3 

Re= 150000, 35400 time steps 

NCC-1 .1 .4. no smoothing 

P-referen ce= 1 OD29G 

2nc=0.0, 4th=Q.0O1 , CFL=1 .0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



2.3100&O4 
2 1005 EM 
2 014ir-£M 
1 U023L-U4 


I 


1 7103C-EM 
1.5S01&W 
1.4O0O&OKJ 
1 .Z530&M 
1 1010PM 
3 490«&O5 
7.9?S3t-0£ 

4 . 9331&05 

3412O&05 



Power spectrum and two-point correlations: 


PS D el ax isl velocity w, P RNS wd h RC P = 0,3 
Pipe llcw, Fe = 150,000 
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Comparison with RANS and experimental data: 


Pipe flow, Re = 150,000 . ... _ . 

^ O PRNS with RCP = 0.3 
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A.3 PRNS with RCP = 0.34 


History of velocity components at probe 1 and probe 14: 


Pipe flow Re=1 50k, at probe 1 
PRNS, RCP=0.34, linear model 



5 


27 r 

26 r 


Pipe flow Re=1 50k , at probe 1 
PRNS, RCP=0.34, linear modeil 
2nd =0.0. 4th =0.001 . cfl=1.0. dt=2.5fr5 
Convergence order=3, Max iter=12Q 


w 



18 0 


i I i i i I i i • i I i i ' i I i i ; I 

1000 0 20000 30000 40000 00000 


Time steps 


Pipe flow Re=15Qk, at; probe 14 
PRNS, RCP=0.34, linear modeil 



Pipe flow Re=1 50k, at probe 1 4 
PRNS, RCP=0.34, linear modeil 
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History of subscale turbulent kinetic energy at probes 1 and 14. 


15r Pipe flew Re= 150k, at probe 1 
14 L PRNS, RCP=Q.34 linear model 

■ 2nd=Q.O. 4th=0.001 p cfl=1 .0. dt=2.5e-5 
13 " Converged oe ord er=3 r Max iter= 1 20 

1.2 j- 

1.1 - * 

i r 


0.9 - 

o.0 V 

0.7 - 
0.6 r 
0.5 ~ 


0.4 t 
0 3 j- 
0.2 r 

0.1 p 

oLl 


.1 . 1,1 I llllLt Ini 1 L . Lt. 


unit i ii^ j 


10000 20000 30000 40000 00000 



Time steps 


Time steps 


Contours of w, u, v, pg, k, /u + ju T at the time step 46,400: 


PR NS, linear model. Rep *0.34 

Re= 150,000, at 46400 time steps 

NCC-t .1 A, no smoothing 

P-referen ce= 1 OD29Q 

2nc=0.0, 4th=D.001, CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



■ 

i 


2.2S70E+O1 

s.uebe+qi 

1 gEOfaETi-Dl 
1 .d274£+0l 
l 6442E«-01 
1 49l0£*0l 
1.3370&-01 
1 1W6&-01 
1 03146*01 
U 70M£*DO 
7.2SOS&-00 
S 710^E1»DU 
4 1007&-OO 
2QM0E+OP 


PR NS, linear model. Rep *0.34 

Re= 150,000, at 46400 time steps 

NCO-t .1 A, no smoothing 

P-referen ce= 1 OD29Q 

2nc=0.0, 4th=Q.001, CFL“1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



I 


u 

2.0633E+OD 
2.2909E+nn 
i f3irasr»nn 
1.540lE*OO 
\ 1737E+00 
B.012TE41 
4 2886 E-Cl 
5.&S43E.02 
.3 1S3T&01 


■H UU30L-01 
-1 .Q 008 E+QQ 


U -l 4332 L*QU 
-1 .0d00e+DD 
.2.178ge+PD 
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PR NS, linear model. Rep =0.34 

Re= 150,000, at 46400 time steps 

NCC-1 .1 A, no smoothing 

P-referen ce= 1 0D29G 

2nc=0 .0 , 4th=Q.0Q1 „ CFL=1.0, dt=2 . 5e-5 

Convergence 3 order or Max iter=120 



Y 

1.001S&-OD 
i .5771 t+on 
i rursnon 
9 . 40 - 10 &O 1 
fj 3370&01 
3 1934&01 
4 . 70 O 0 &O 3 
- 3 .Q 9 Q 8&01 
s rflurni 
-S. 3 M 0 &O 1 
-1 2 S 31 E+ 0 Q 

I ” " . 5075 E *00 
*1 .002QE+QD 
. 2 . 1954 &-DD 


PR NS, linear model, Rep =0.34 

Re= 150,000, at 46400 time steps 

NCC-1 .1 A, no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=Q.0O1 „ GFL=1 .0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



Pa 

0.0347&-O0 
5.9is^E+nn 
1 s ig^nEiton 

3.Q3eQE*DQ 
2.3107E+00 
.5594E+DQ 
0 0 OO 2 E -0 1 
! hi U hi IJ LA 1 
-S .50056-01 
-1 27UUL»UU 
1 .1 .9073&-0Q 
I -2.7H37EM1D 




PRNS, linear model. Rep =0.34 

Re= 150,000, at 46400 time steps 

NCC-1 1 4, no smoothing 

P-referen ce= 1 0D290 

2nc=Q.O, 4th=D.0O1 „ CFL=1 .0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



luJt^ 

.oooob-oo 

.0Q77E+DD 
01WE*OD 
OlE*QQ 
jfi :■ r&tii, hjij 
2,0S85E*QQ 
l.MttE+00 
1 .0050€+OD 
i «isr*nn 
1 .2dS2&-DQ 
1 ,Q700E*QQ 
e.M02E-ai 
0.0201 SOI 
s.aaaa&oi 



PRNS, linear model. Rep =0.34 

Re= 150,000, at 46400 time steps 

NCC-1 1 4, no smoothing 

P-referen ce= 1 0D290 

2no=0.0, 4th=D.0O1 „ CFL=1 .0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



I 

I 

■ 


mu 

5.0<110&(M 
5.53iaS6M 
£ 12WE4H 
4 . 7 Q 99 &M 
4 23 B &04 
3 UIJMIJ L-0-S 
3 4703 &M 
3.QS77&04 
2 0S71&CM 
2.MH&M 
1 0300 &(M 
1 .425SE-(H 
1.01406-04 
5.Q440&Q5 



A.4 PRNS with RCP = 0.4 

History of velocity components at probe 1 and probe 14: 


Pipe Flav, Re= 150000, at Probe 1 
PRNS, RCP=Q 38 

2nd=Q.O, 4th= 0.001, CFL=1, dt=2 5e-5 
Co iwergen ce order= 3, or M ax. iter= 1 20 



10000 15000 20000 25000 30000 

Time Steps 
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Time Steps 



History of subscale turbulent kinetic energy at probes 1 and 14: 
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Contours of w, u, v, pg, k, /u + ju T at the time step 29,000: 


PRNS, Rep =0.38 

Re=1 50000 , 30000 time steps 

NCC-1 .1 .4, no smoothing 

P-referen ce= 1 0D29G 

2nc=0 .0 , 4th=Q.0O1 p CFL=1.0, dt=2 . 5e-5 

Convergence 3 order or Max iter=120 



I 


2 11866*01 
1 .97376+01 
1 03006N01 
1.00796+01 
1 WS0E-01 
1 JU22DU1 
1.25936*01 
1 11046+01 
0 73576*00 
0.9O7OE+00 
5?7?:?E+0O 
5 4430L*0U 
4.O2O0E+00 
2,5922 6+00 


PRNS, Rep =0.36 

Re= 150000, 30000 time steps 

NCC-1 .1 .4. no smoothing 

P-referen ce= 1 0D29G 

2nc=0.0, 4th=D.Q01, CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



0.2250601 
5.30S0&Q1 
4 30526=01 
3 4049 601 


2 54456m 


I 


1 .83426-01 
7 .03006-02 
■2.1055&02 
1 1300101 
-2.05736-01 
“1770E-1 1 
-3 030UL-O1 
-4.0103 601 
•5.73076-01 



PRNS, Rep =0.36 

Re= 150000, 30000 time steps 

NCC-1 .1 4. no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=D.001 „ CFL=1 .0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



i 


1 .22346*00 
1 .04576+00 

0 070£f-Q1 
9.9022&01 
5 12486-01 
3 3475&01 
1 .57026-01 

-2 ■07176-02 

1 9MSF01 
-3 .78106-01 
-5.5392&01 
-7 .3 10S6-O1 
-0.0938601 

■ 1 09716+00 


PRNS, Rep =0.36 

Re=150000, 30000 time steps 

NCC-1 1 A, no smoothing 

P-referen ce= 1 00290 

2no=0.0, 4th=D.001 „ CFL=1 .0, dt=2.5e-5 

Convergence 3 order or Max tter=120 



■ 


p'j 

3.05206*00 
3.1929 6+QQ 
2 0337 ft *00 
2 474&LMJU 
2 11536*00 
1.75016*00 
1.39896*00 
1 .03776+00 
B 7055601 
3.1937601 
-3 .98 17 602 
-3.93UUL-01 
-7.5010&O1 
1.11746+00 



PRNS, Rep =0.36 

Re= 150000, 30000 time steps 

NCC-1 .1 4. no smoothing 

P-referen ce= 1 00290 

2nc=0.0, 4th=0.QO1, CF 1=1.0. dt=2 5e-5 

Convergence 3 order or Max iter-120 



I 

I 

i 


iuitv 

3.50706+00 
3 30356+00 
3.09016*00 
2.00406*00 
2 09LHL' UU 
2.40016*00 
2 20176+00 
2 07746*00 
1 .07306*00 
1 05076*00 
1 -10436+00 
1 20006+00 
1 05506+00 
0.51306-01 


PRNS, Rep =0.36 

Re= 150000, 30000 time steps 

NCC-1 1 4. no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=D.001 , CFL=1 .0, dt=2.5e-5 

Convergence 3 order or Max tter=120 



I 

I 

■ 


2 40056=03 
2.2072603 

2 1130603 
1 WUSL-Q3 
1 70716-03 
1 &93UL-(f3 
1.4204603 
1.24716-03 
1 0737 T 03 

3 00406-04 
7 2705604 
5.53706-04 
3. 0035 604 
2.070O&O4 
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A.5 PRNS with RCP = 0.5 


History of velocity components at probe 1 and probe 14: 


Pipe Flow, Re= 150000, 3k Probe 1 





Pipe Flow. Re= 15OOO0, at Probe 14 
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History of subscale turbulent kinetic energy at probes 1 and 14. 


Pipe Flow, Re= 150000, £ Probe 1 




Contours of w, u , v, pg, k, ju + ju T at the time step 25,000: 


PRNS h Rep =0.5 

Re- 150000, 25000 time steps 

NCO-11 4. no smoothing 

P-referen ce- 1 0D290 

2no-Q.O, 4th-0.001 „ CFL-1 .0, dt-2 5e-5 

Convergence 3 order or Max tter=120 



I 


1.015OE+Q1 
1.09316*01 
i S713H01 

13275E*0l 


0 

7. 

5. 



1 


4.74326*00 

3.02436*00 

2.30556+00 


PRNS, Rep =0.5 

Re- 150000, 25000 time steps 

NCO-1 1 A. no smoothing 

P-referen ce- 1 0D290 

2nc-Q.O, 4th=D.001 „ CFL-1 .0, dt-2 5e-5 

Convergence 3 order or Max tter=120 



■ 


5. 2 130601 

4.3005601 

3 52oair-m 

26735601 
l 02TQ6P1 
9.0040602 
1.3390602 
7 1253602 
1 5£90F (11 
*2 4055601 
-3 252Q601 
4.0905601 
4.9450 601 
'5 7915601 



PRNS, Rep =0.5 

Re-15O0OG, 25000 time steps 

NCO-11 4. no smoothing 

P-referen ce- 1 0D290 

2nc-Q.O, 4th-0.001 „ CFL-1 .0, dt-2 5e-5 

Convergence 3 order or Max tter=120 



I 


K 

1 .20236*00 
1 .09546+00 
h nr; 'fir si 
7.2159601 
5 34 "Ol.m 
3.4701601 
1.0092601 
•2.5970602 
■2 1:0017(11 
*3.9975601 
-5.0004601 
-7 7353 601 
-9.0042601 
• 1.14736+00 


PRNS, Rep =0.5 

Re-150000, 25000 time steps 

MCC-1 .1 4. no smoothing 

P-referen ce- 1 0D290 

2nc-Q.O, 4th-D.001, CFL-1 .0, dt-2 5e-5 

Convergence 3 order or Max Eter=120 



i 


3 7944 E+Q0 
347096+00 
3 1593 El* DO 
2. 04 10 £*00 
2 52426*00 

2 200 7 L 1 00 
1 .0092 E+DO 
1.57106*00 
1 2£4 16* 00 
9 305 7 60 1 
0 1904601 

3 0151601 
-1 .0025602 
■3.3350601 
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U, V 


PRNS, Rep =0.5 

Re=1 50000 , 25000 time steps 

NCO-1 .1 .4, no smoothing 

P-referen ce= 1 00290 

2nc=0 .0 , 4th=O.OG1 , C F L= ' 1 . 0, dt=2 . 5e-5 

Convergence 3 order or Max iter=120 



I 


lurtv, 

2.soaaetQD 

2.31S4E+DD 
2 lanartDn 
l.M93EtQD 
l JGlSfii'Qu 
1 5?M£*£lU 
1.30236*00 
1 .2Q77€+Q0 
i a:3ir*0Q 
0.3M9&O1 

a.ssese-oj 

4 FJ323&-01 
2 .3432641 
i .aaga&oi 


PRNS, Rep =0.5 

Re=1 50000 , 25000 time steps 

NCO-1 .1 .4, no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=Q.0O1 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



trru 

2 7144&Q3 

?M17rM 

2.14M&03 


I 


l 049OE-O3 
1 7577fr&3 
1.5883643 
i .37505 m 
1 103(3f m 

8 QQ02E-C- 
3.0953&D4 
4 1021SO4 
2.2090&SM 



A.6 PRNS with RCP = 0.8 


History of velocity components at probe 1 and probe 14: 


Pipe Flow. Re= 150000, at Probe 1 
2 PRNS, RCP=0.S 

2nd=D.0, 4th=D.O01, CFL=1, dt=2 5fr 5 
Co rwergan ca order= 3, or M ax itsr= 1 20 


u 



-i r 


-2 ■ 

1 5001 10001 15001 20001 

Time Steps 


Pipe Flow. Re= 150000, at Probe 1 
25 r— PRNS, RCP=0.8 

2nd=Q.0, 4th= Q .00 1 , CFL= 1 . dl=2 5e-S 
^ Co iwergen ce order= 3, or Max. iter= 1 20 



16 r 


0 5000 10000 15000 20000 25000 

Time Steps 


Pipe Flow, Re= 150000, d Probe 14 
2 r- PRNS, RCP=0.8 

2nd=0.0, 4th= Q .00 1 , CFL= 1 . dt=2. 5e-5 
Co iwergen ce order= 3, cr Max. iter= 1 20 



Time Steps 


Time Steps 
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History of subscale turbulent kinetic energy at probes 1 and 14. 



Pipe Flow. Re-ISQODO, at Probe 14 
7 r PRIMS, RCP=0 S 

2nd=Q.O, 4th= Q .00 1 , CFL= 1 . dt=2 5e-5 
Co iwergen ce order =3, or M ax iter= 1 20 



Time Steps 


Contours of w, u, v, pg, k, // + ju T at the time step 25,000: 


PRNS, Rep =0.8 

Re= 150000, 25000 time steps 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 0D29Q 

2nc=0.0, 4th=D.001 p CFL=1.0, dt=2.5e-5 

Convergence 3 order or Max iter=120 



I 


i. 0 saae+oi 

1.45556*01 
1 351 1E+05 
1 .24076+01 
i njaf-nji 
1 03796+01 
9. 3S53&-00 
0.29146+00 
7 2474E+00 
0 .20356+00 
5 . 1 5956*00 
4 11506+00 
3.07106+OQ 
2. 02776+00 


PRNS, Rcp=0.8 

Re= 150000, 25000 time steps 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 OD29G 

2nc=0.0, 4th=Q.001, CFL“1.0. dt=2 5e-5 

Convergence 3 order or Max tter=120 



I 


1 .0709 E+00 
0.0005601 
0 5435601 
442056-01 
2 29756-PI 
1 .74506-02 
-1 0405601 
-4 07156-01 
JJ 13456-01 


-0 ,3 1 75 L-0 1 
-1 .04406+00 


1 -1 .25336+00 
-1 4000 6+0D 
1 00096+00 



PRNS, Rcp=0.8 

Re= 150000, 25000 time steps 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 0D29G 

2nc=0.0, 4th=D.001 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 


I 


X 

1.07016+00 
1 40 12 E+00 
1 1 234 El ► 00 
0.455 7601 
5 0773601 
2.69906-01 
1 20076-02 
•2 .0577601 
5 43006-01 



PRNS, Rcp=0.8 

Re= 150000, 25000 time steps 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 0D29G 

2nc=0.0, 4th=D.001 p CFL-1.0, dt=2.5e-5 

Convergence 3 order or Max tter=120 



I 


pg 

4.50306+00 
4.15226+00 
3 72156+00 
3.29076+00 
2 05006*00 
2.42926+00 
1 .09046+00 
1 .50776+00 

1 1300E+OO 
7.00156-01 

2 7539601 
-1 .5537601 
-5.0013601 

• 1.01096+00 
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PRNS, Rcp=O.S 

Re=1 50000 , 25000 time steps 

NCO-1 .1 .4, no smoothing 

P-referen ce= 1 00290 

2nc=0 .0 , 4th=Q.0G1 , C F L= ' 1 . 0, dt=2 . 5e-5 

Convergence 3 order or Max iter=120 



I 


luftn, 

3.5O70M 

3.3035e+nn 

3 uggie*0D 
2 .0940**00 
3 syrwi i uii 
2.40016*00 
2 20176+00 
2.0774E+Q0 

i 07306*00 
1.00076+00 
1 .40436+00 
1 20OQIL+OU 

i .ossae+oo 

0.513O&01 


PRNS, Rcp=0.£ 

Re=1 50000 , 25000 time steps 
NCO-1 .1 .4, no smoothing 
P-referen ce= 1 00290 
2nc=0.0, 4th=D.001 p CFL=1.0, dt=2 5e-5 
Convergence 3 order or Max iter=120 



mu 

3.3201603 
3.09326 03 
2 05046-03 
2 .0235603 


3 300F3E.U3 
2.15376-03 
1. SI 09603 
1.WD&Q3 
1 44016 03 
I 21421.-03 
0 7930 604 
7 44406-04 
5. 0900 60a 
2 . 7 - 17360-1 



A.7 PRNS with RCP = 1.0 


History of velocity components at probe 1 and probe 14: 


Pipe Flow. Re= 150000, at Probe 1 
2 r- PRNS, RCP=1.G 

2nd=0.0, 4th=0.001, CFL=1, dl=2 5e-S 
Co nvergen ce order= 3, or M ax iter= 1 20 



-1 . 


1 i i i i I i i i i I i i i i I i i i i I i i i i I i 

0 10000 20000 30000 40000 50000 

Time Steps 


Pipe Flow. Re= 150000, d Probe 1 
25 tr PRNS, RCP=1.G 

2nd=Q.0, 4th= 0.001, CFL=1, dl=2 5e-S 
24 f Co iwergen ce order= 3, or M ax <ter= 1 20 



0 10000 20000 30000 40000 50000 

Time Steps 



Pipe Flow. Re= 150000, at Probe 14 
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History of subscale turbulent kinetic energy at probes 1 and 14. 



Time Steps 


Pipe Flow. Re= 150000, at Probe 14 



Contours of w, u, v, pg, k, // + ju T at the time step 54,000: 


PRNS, Rcp=1.0 

Re= 150000, 54000 time steps 

NCO-t .1 .4. no smoothing 

P-referen ce= 1 0D29Q 

2nc=0.0, 4th=D.001 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



I 


1.3347M1 
1.23W&-01 
i ugie+os 

1 Q013E+0J 
a 73556 >pu 
UU570LHJU 
7.37976+00 
7 1 □ 1 B E+DQ 
15 2230E+DP 
5 .34906+00 
4.4901 6+00 
3.5302E+O0 
2.71236+00 

1 9344 E+0Q 


PRMS, Rcp=1.0 

Re= 150000, 54000 time steps 

NCO-1 .1 .4. no smoothing 

P-referen ce= 1 OD29G 

2nc=0.0, 4th=Q.001, CFL“1.0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



I 


9.SS99&01 
5.0512&O1 
4BUH&01 
3 7545E-01 
2 TCn^&Ol 
1 4570L-01 
9.D9S0&02 
4.30056,02 
1 4073&O1 
-2.53SS&01 
-3.59396-01 
4.9322&01 
-5.98flfi6.01 
.0.7299&01 




PRNS,Rcp=1.0 

Re= 150000, 54000 time steps 

IMCO-1 .1 .4. no smoothing 

P-referen ce= 1 0D29G 

2nc=HO, 4th=0.001 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max iter=120 



9.03196.01 
0.3127&O1 
S 793£f-01 
5 27436-01 
3 T552&01 


2.2300&01 
7. 1003 6-02 
■S .02335-02 
2 321Sf 01 
-3.0407&O1 
-5 35086-0-1 


I 


-hi 07906-0 I 
-0 3902&01 
■9.9173&01 


PRNS,Rcp=1.0 

Re= 150000, 54000 time steps 

NCO-t .1 .4. no smoothing 

P-referen ce= 1 0D29G 

2nc=0.0, 4th=0.001 p CFL-1.0. dt=2.5e-5 

Convergence 3 order or Max tter=120 



I 

i 

i 


4 1574 5+00 
3.91016+00 
47rg non 
3 1307E+00 
2 7005E+OQ 
2 44026+00 
2 10406+00 
1 70196+00 
1 4105E+OO 
1.07736*00 
7.35006-01 
3.9200i-Q1 
5.00276.02 
.2.9100501 
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PRNS, Rcp=1.Q 

Re=1 50000 , 54000 time steps 

IMCO-1 .1 .4, no smoothing 

P-referen ce= 1 00290 

2nc=0 .0 , 4th=O.OG1 , C F L= ' 1 . 0, dt=2 . 5e-5 

Convergence 3 order or Max iter=120 



I 


lurt^ 

2.00006+00 

i .gggye+DD 
i nn’soDn 
1 .71 126+O0 

l hl^MI i Nil 
1 5 1 IJhi L 1 UU 
1 .42346+00 
1.33016*00 
i 22S8 e+oo 

1 1330E+QO 
1.0»T3E*00 
y^iQsE-ai 
0 4470&O1 
7 40406-01 


PRNS, Rcp-I.Q 

Re=1 50000 , 54000 time steps 

MCC-1 .1 .4, no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=Q.0O1 p CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max tter=120 



S .95066.03 
5.63S06-Q0 
S 1130IT-O1 
4,4914£.03 


■I 4 20091M13 

■ 3 6404&O3 
3 J239E-03 

■ 3 00 ICE 03 
m 3 0T80E-03 
m 2,1535£-Q3 

1 73406-03 

I I 3114L-03 
0.09006.04 
4,g001&O4 



A.8 Steady RANS simulation 


Iteration history of velocity components at probe 1 and probe 14: 


Pipe Flow, Re=15OD0O, at Probe 1 



Pipe Flow, Re=15OD0O, at Probe 1 



Pipe Flow, Re=15OD0O, at Probe 14 
Steady RANS Model with WF 
2nd = -0.01 ,4tti = D.05,CFL = 3 


Pipe Flow, Re=15OD0O, at Probe 14 
Steady RANS Model with WF 
2nd =-O.01, 4th =0.O5, CFL = 3 


o - 


-003 


U 

v 


f 


1 ■ ■ 1 * 1 * 1 1 1 1 1 1 1 * • • 1 1 I 

□ 100000 200000 300000 400000 


w 



° 0 100000 200000 300000 400000 


Pseudo-time iterations 


Pseudo-time iterations 
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Iteration history of turbulent kinetic energy at probes 1 and 14. 




Contours of w, u, v, pg, k, // + ju T at the iteration number 400,000: 


Steady RANS with WF 
Re=1500C0, iterations; 400000 
NCO-1 .1 4. no smoothing 
P-referen ce- 1 0D290 
2nd— 0 01, 4th =0 05, CFL-3.0. 



i.Msse*oi 

i . 24506+01 
1 1BB1£*D1 
1. 00048*01 

1 00F5lgi.pi 

9 .2TD2E*00 
0.47326+00 
7.07026+00 
B 07928+00 
0 00236* UU 
5 2053 E+QQ 

1 4 40036+00 
3.09136+00 
2 99446+00 


Steady RANS with WF 
Re-150000, iterations: 400000 
NCO-1 1 4. no smoothing 
P-referen ce- 1 00290 
2nd— 0.01, 4th-0 05, CFL-3 0. 



I 


0.0010600 
7 .0920600 
B 1239F-0B 
5,1540605 
4 10506-00 
3.2100600 

2 2470600 
1 .2700600 

3 0900 6 07 
-0.0UOl6.07 
-1 .0290609 
-2 .5900600 
.3,50706.00 
.4.5300600 



Steady RANS with WF 
Re-150000, iterations: 400000 
NCO-1 .1 4. no smoothing 
P-referen ce- 1 00290 
2nd— 0.01, 4th-0 05, CFL-3 0. 



I 


V 

1 11006.02 
94000603 
7 7450603 
0.0240603 
4 30376-03 
2.50276-03 
0.0173604 
■0.5023604 
-2 50O2F 03 
4.30126-03 
-0.0221 603 
-7 7431 603 
-9 4041603 
1.1105602 


Steady RANS with WF 
Re-150000, iterations: 400000 
NCO-1 1 4. no smoothing 
P-referen ce- 1 0D290 
2nd— 0.01, 4th-D 05, CFL-3 0. 



i 


1.2902EtQd 
1 19096+00 
1 05206+00 
9 

0 2000601 
7.0777601 
5.0940601 
4 7110601 
3 5290601 
2.3401E-O1 
1.1032601 
-1 .9000 603 
.1 .2020601 
.2,3956601 
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turb* u, v, w 


Steady RANSwithWF 
Re= 150000, iterations; 400000 
NCO-1 .1 .4. no smoothing 
P-referen ce= 1 OD29Q 
2nd— 0.01, 4th =0 05, CFL=3 0. 



I 

■ 


lurtv, 

1 .63306*00 
1 . 50226*00 
i 4gogrnDO 
1 4190 6*0U 
1 34031 i. uu 
1.2?Q9E*DQ 
1.20606*00 
1 13436*00 
1 053SE+00 
3.91006-01 
9 2027601 
0 40346-O1 
7 7701&O1 
7. 00236-01 


Steady RANS with WF 
Re=1 50000 , iterations; 400000 
MCG-1 .1 .4. no smoothing 
P-referen ce= 1 00290 
2nd— 0.01, 4th =0.05, CFL=3 O. 



1 


mu 

5,03306.03 
4.0701603 
4 3193603 
3 9024603 




3 KB0&03 
3.24076-03 
2.0919603 
2.5350603 
2 1701&03 
1 0213L-03 
1.4044603 


I 


1 1070603 

7.5072604 

3.9300604 




Radial profiles of w, u, v, k and ju + p T : 


15 - 


Steady RAMS with WF 
R 6=150000, iterations: 400000 
NCC-1 .1.4, no smoothing 
P-reference= 1 00290 
2nd=-0.01, 4£h=0.O5, CFL=3.D. 



u 

V 


w 


0.O4 r- 


0.03 - 
0j 0£ r 


Steady RAMS with WF 
Re= 150000. iterations: 400000 
NCC-1 .1.4, no smoothing 
P- references 1 00290 
2nd=-0.Q1. 4th=0.O5, CFL=3.D. 


u 

v 


-0.01 - 


0 001 0.02 005 0.04 0.05 0 06 0.07 

Y 


0.05 0.04 

Y 


0.05 0.06 0.07 


2.5 r 


Steady RAMS with WF 
Re=15OQ00, iterations: 400000 
NCC-1 .1.4, no smoothing 
P- references 100290 
2nd=^D.01 , 4th=0.O5, CFL=3.D. 


lutb^ 






0 


0 


_L_I_I ■ ■ I ■ ■ I » ' ■ I ■ • 1 ■ I ■ ■ ' I I ■ I I 

0 01 0.02 0 05 0.04 0.05 0 06 0.07 

Y 
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Appendix B: Effects of the Resolution Control Parameter on the PRNS 
Solutions (Re= 15,000) 

B.l PRNS with RCP = 0.0 

History of velocity components at probe 1 and probe 14: 


Pipe Flav. Re= 15000, at Probe 1 

DNS, 2nd=D.O, 4lh=D.0Q1, CFL=1, d=4.eh4 

Coiwergan cs order= 3, or M ax it»r= 1 20 





Time Steps 


Pipe Flow, Re= 15000, at Probe 1 
4 p- DNS, 2nd= D O, 40i= D.0D1 , GFL= 1 , d=4.e-4 
Cojwergan ce order= 3, or M ax iter= 1 20 



1.5 - 


1 - 



I 1 i i i I i i I i i i ■ I 

2500 5000 7500 10000 

Time Steps 


Pipe Flav. Re= 15000, at Prnbe 14 



Pipe Flow. Re= 15000, at Prnbe 14 
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Contours of w, u, v, pg at the time step 10,000. 


ONS. Re=15000. 10000 time steps 
NCC-1 .1 A, no smoothing 
P-referen ce= 1 00290 
2nc=0.0, 4th=0.001, CFL=1.0, clt=4e^4 
Convergence 3 order or Max tter=120 



I 2.52(HE*00 
2.34146+00 
2 i024F»Dn 
1.90341i»UU 

■ 1 BtMtt, HJIJ 
1.0256LMJU 
1 4405 E+00 
1 .20756+00 
i ana^rnna 
y uaiOL-tfi 
7.3001 601 
5.5103601 
3 7205601 
1.9309601 


ONS. Re=1 5000, 10000 time steps 
NCC-1 .1 A, no smoothing 
P-referen ce= 1 00290 
2nc=0.O, 4th=0.001, CFL=1.0, clt=4e^4 
Convergence 3 order or Max tter=120 



I 

i 


2 7525603 
2.4009603 
2 0013603 
1.7157603 
1 3101603 
1 .0245603 

0 7905 604 
3.3325604 

1 2300F flS 
-3 5737604 
-7.0357604 
-1 IP32L-03 
.1 3940603 
.1 7404603 




DNS, Re=1 5000, 10000 time steps 
NCC-1 .1 A, no smoothing 
P-referen ce= 1 00290 
2nc-0.0, 4th=0.001, CFL=1.0, dt=4fr4 
Convergence 3 order or Max tter=120 



1 .00056 03 
9.3343604 
0 55:i-rn-: 
4.9931604 
3 3225604 
1.0519604 
-1 .0700 E’00 
1. 0993 E 04 
■3 3599 FEW 
-5.0305604 
-0 7011604 

1 -0.3717604 
.1 .0042603 
■ 1.1713603 


ONS. Re=1 5000, 10000 time steps 
NCC-1 .1 A, no smoothing 
P-referen ce= 1 00290 
2nc=0.0. 4th=0.001 , GFL=1 .0, dt=4fr4 
Convergence 3 order or Max tter=120 



i 


2.0200604 
2.4939604 
2 1591604 
1 .0243 604 
l 488S&04 
1 1540 6EM 
0.1301605 
4.9499605 
1 5010FOS 
-1 .0407605 
-5 1349605 
4.5432605 
.1 1091604 
.1 5240 604 



B.2 PRNS with RCP = 0.2 

History of velocity components at probe 1 and probe 14: 


Pipe Flow, Re= 15000, at Probe 1 
0 2 i- PRNS, RCP=0 2 

2nd=Q.0, 4th= 0 .00 1 , CFL= 1 . dt=4. e-4 
Co iwergen ce order= 3, cr Max. iter= 1 20 



QJ05 r 


> 



0 - 


-OJ05 - 



-o.i r 

■0.15 « 

-0.2 — 1 * 1 1 1 1 1 1 1 ! 1 1 

0 5000 1 0000 

Time Steps 


Pipe Flow, Re= 15000, at Probe 1 
5 r- PRNS, RCP=0.2 

2nd=0.0, 4th= 0.001, CFL=1, dt=4. a-4 
^ Co iwergen ce order= 3, or M ax <ter= 1 20 




3 5 r 
3 - 

2.5 r 

2 - 

1.5 r 

1 r 

0.5 r 

Q'— 



Time Steps 
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Pipe Flow. Re=15QO0, at Probe 14 
2 r PRNS, RCP=0 2 

2nd=0.Q, 4th= 0.001, CFL=1. dt=4.a-4 
Co iwergen ce order =3, or M ax iter= 1 20 



J ! i i I i i I | i J 

0 5000 10000 15000 

Time Steps 


History of subscale turbulent kinetic energy at probes 1 and 14: 


Pipe Flow. Re= 15000, at Pmbe 1 
PRNS, RCP=0.2 

2nd=Q.0, 4th= 0.001, CFL=1. dt=4.a-4 
Co iwergen ce order= 3, or M ax iter= 1 20 


5000 1 0000 

Time Steps 


002 


OjOI 


Pipe Flow. Re= 150O0, at Pmbe 14 
PRNS, RCP=0.2 

2nd=Q.0, 4th= 0.001, CFL=1, dt=4.a-4 
Co iwergen ce order= 3, cr Max. iter= 1 20 


50 00 1 0000 

Time Steps 
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Contours of w, u, v, pg, k, /u + ju T at the time step 15,000: 


PRNS, Rep =0.2 

Re= 15000, 15000 time steps 

NCC-1 .1 A , no smoothing 

P-referen ce= 1 00290 

2nc=G.O, 4th=0.(J01 . CFL=1 .0, dt=4e-4 

Convergence 3 order or Max iter=120 



n 

I 2.98S2&-00 
2.75796+09 
2 M7S£ »DD 
2.3372E+9Q 
.■ nm mu 
i,3i03H*uu 
i Toaie+QQ 
1 49506+99 
1 20WE+DD 

i.nsi£*un 

8.W7SE-01 
4.W4UL-01 

4 «oa&oi 

2.0071601 


PRNS, Rep =0.2 

Re= 15000, 15000 time steps 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=D.Q01, CFL=1.0, dt=4e-4 

Convergence 3 order or Max £ter=120 



I 

i 


T.saoa&tM 

9 4001 604 
5 M14T £M 
4.2l00[i.EM 
3 D919E.EM 
i yy l C-tM 

0 4230605 
■2 .0200605 

1 4971604 
-2.S3l3ti.tM 
-a.assa&tM 
-4 7014&IM 
- 5 . 990160 * 

•7.9099604 



PRNS, Rep =0.2 

Re= 15000, 15000 time steps 

NCC-1 .1 .4. no smoothing 

P-referen ce= 1 0D290 

2nc=0.0, 4th=D.001, CFL=1.0, dt=4e-4 

Convergence 3 order or Max iter=120 



■ 3.2029604 

■ 2.75396 04 

■ 2 3952604 

Q 1 0503604 

1 flUHlf -EM 

u ijuuuL-es 
S.1109&O5 
9.33076.09 
.3 0500 F OS 
-U .3347605 
-1 2019&O4 
-1 73926fM 
-2 1700604 
•2 9270 E- 04 


PRNS, Rep =0.2 

Re= 15000, 15000 time steps 

NCC-1 .1 .4. no smoothing 

P-referen ce= 1 0D29Q 

2nc=0.0, 4th=Q.001, CFL=1.0, dt=4e-4 

Convergence 3 order or Max iter=120 



I 

i 


P9 

6 3052605 
9.9719605 
9 77076-05 
U JUWL-tfS 
0 182 1605 
7 .0903603 
7 .09506 05 
7.3134605 

7 aiaii; ns 
0 .7250605 
9 4320605 
0 1333605 
5.9401 60S 
5.5520605 



PRNS> Rep =0.2 

Re= 15000, 15000 time steps 

NCC-1 .1 4. no smoothing 

P-referen oe= 1 00290 

2nc=0.0, 4th=D.0O1 „ CFL=1 .0, dt=4e-4 

Convergence 3 order or Max iter=120 



I 

■ 


1 .00946 09 
1 .5514609 
1 43356 09 
1.31556 09 
1 1970603 
1.979011-03 
9 0197619 
0.49726-10 
7 2570610 
9 .070361 D 
4 0909619 
3.71 33 L- 1 U 
2.5399619 
1.39046 1 9 


PRNS> Rep =0.2 
Re=15000, 15000time steps 
NCC-1 1 4. no smoothing 
P-referen ce= 1 00290 
2no=0.0, 4th=D.0O1 „ CFL=1 .0, dt=4e-4 
Convergence 3 order or Max iter=120 



I 

■ 


1.0079 60S 
1.9979605 
I 0979 T 75 
1 M7UL-Q& 
1 .0070605 
1.00 79 60S 

1.0079605 
1 .9979605 
1 007017 ns 
1.0070605 

1 .0079605 
1 .0079605 
1.9979 60S 
1 .9979605 
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B.3 PRNS with RCP = 0.3 


History of velocity components at probe 1 and probe 14: 




Pipe flew Re=15k, at probe 14 
PRNS, linear model, RCP=0.30 
2nd=Q.Q„ 4th=O.QOOQ1 , cfi=1.0, dt=4.Qe-4 
Convergence order=3. Max rter=120 

li 


Pipe flew Re=15k, at probe 1 4 
PRNS, linear model, RCP=0.30 
2nd=Q.Q. 4th=O.QOOQ1 , cfl=1.0, dt=4.Qe-4 
Convergence order=3. Max rter=120 
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History of subscale turbulent kinetic energy at probes 1 and 14. 




Contours of w, u, v, pg, k, ju + p T at the time step 15,000: 


PR MS, linear model, Rop-0.30 
Re= 15,000, at the 15,000 time step 
MCC-1 1 4. no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 0, ctt=4 Qe-4 
Convergence 3 order or Max iter=120 



I 


i 


¥¥ 

1.7&1B&-09 
1 .593 16+99 
1 40446*99 
1.34S7E+99 
1 22?iE+09 

1 1U04L*UU 
9 .0971601 
0 71046-01 
7.5239601 
0.3390601 

5.1591601 

3.9033601 

2 7785 601 
1.5090601 


PR MS, linear model, Rcp=Q.3Q 
Re= 15,000, at the 15,000 time step 
MCC-1 .1 A. no smoothing 
P-referen ce= 1 0D290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Qe-4 
Convergence 3 order or Max iter=120 



■ 4.5S7Q&02 

■ 3.017B&Q2 

■ 3 0490602 
2 2700602 


I 

I 


1 5195602 
Ml 39603 
-2 7737&04 
7 .00006 03 
1 £880r07 
-2.3351602 
-3.1042602 
-3 .0734602 
-4 8425602 
-5.4118602 



PR MS, linear model, Rcp=0.3Q 
Re= 15,000, at the 15,000 time step 
MCC-1 .1 A. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max iter=120 



I 


0.0025602 
0.1W2&02 
E 3255602 
4 5475602 
3 7093607 
2.99996 02 
2.2120602 
1 4342602 
0££g3rO3 
-1 .2249603 
-3 .9972 603 
-1 .8739 6 02 
-2 4574602 
-3.2357602 


PR MS, linear model, Rop=0.30 
Re=15 n OO0, at the 15,000 time step 
MCC-1 .1 A, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 0e-4 
Convergence 3 order or Max iter=120 



I 


po 

1.9301602 
1.700 1&02 
1 4741602 
1 .24226 02 


i 


l □ 107607 
7 70106 03 
5 4019603 
3 1421603 
0 2229 1! £M 
-1 4977603 
-3 0175603 
-0 1374603 
-0 4573603 
• 1.Q777&02 
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PR NS, linear model. Rep *0.30 
Re=15,000, atthe 15,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 0, ctt=4 Oe-4 
Convergence 3 order or Max iter=120 



I 


Lub, 

2.9757&02 

r^w-o: 1 

2.200Q&O2 

3 0« lE-02 
1 .0SO2&O2 
1 &SM602 
1 4G35&02 
1 340I5£ 02 
1 LWr£-U2 
6. 40056-03 
0 3S87&H3 

4 33Q9&03 
1.1322 & 03 


PR NS, linear model. Rep =0.30 
Re= 15,000, at the 15,000 time step 
NCC-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 0, cft=4 Oe-4 
Convergence 3 order or Max tter=120 



mu 

2.2721&IM 


I 


2.1232&CM 
1 0743 IT- £M 
1 .U2W L-tK 
1.OT064M 
1.S3TQ&M 
1 3707&M 
1 .2200&W 
1 9000604 

a .01091-05 

7 93Q0E-C5 
hi 'J41UL-05 
4.8S20 & m 
3.3S37&9S 



B.4 PRNS with RCP = 0.34 

History of velocity components at probe 1 and probe 14: 




NAS A/TM— 2008-2 15418 


65 


Pipeficftv Re=15k, at probe 14 
PR NS, linear model, RCP=0.34 
2nd =0.0. 4th=G.Q0001 , cfl=1.0, dt=4.0e-4 
Convergence order=3. Max iter=120 

u 


Pipe flew Re=15k, at probe 14 
PR NS, linear model, RCP=0.34 
2nd =0.0. 4th=0.00001 , cfl=1.0, dt=4.0e-4 
Convergence order=3, Max rter=120 
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History of subscale turbulent kinetic energy at probes 1 and 14: 




NAS A/TM— 2008-2 15418 


66 


Contours of w, u, v, pg, k, /u + ju T at the time step 15,000. 


PR NS, linear model. Rep »0.34 
Re=;15,000 p at the 16,000 time step 
NCC-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max iter=120 



I 


1 .2662&-0D 
1.1 7925+00 
1 OBUS+DD 
1 UUtUL'UU 
9 1793601 
0 .3000601 
7.4379601 
0.5072601 
5 60645 01 
4 0257601 
3.9550501 
3 U042I-01 
3.2135&01 
1 3429 501 


PR NS, linear model. Rep =0.34 
Re= 15,000, at the 16,000 time step 
NCC-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



I 


2.07705O5 

2.5039605 

2 lamtr-os 

1 7504605 
i 392U&05 
1.0000605 
0.3508600 
2 .6132600 
1 13+10 00 
-4.0020600 
-8 eooeE-ce 


I 


*1 .2337505 
,1 .0075505 
1 99125.05 



PR NS, linear model, Rop=0.34 
Re= 15,000, at the 16,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max iter=120 



I 


4.9892502 
4.22005-02 
3 4E34F-Q2 
2.0000602 
l 9108602 
1 1512602 
3 03706 03 
-3 .03036-03 
1 1E10F 07 


*1 .3104602 
-2.0058602 

1 -3 4532 L-02 
4.2200502 
4.0000502 


PR NS, linear model. Rep =0.34 
Re=15,000, at the 16,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=Q .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



I 


w 

1.63536 02 
1.3719502 
1 2003 F- 02 
1 U44UL-02 
0 0135&O3 
7.1 7006-03 
5.5437503 
3.9009503 
3 27306.03 


hi yyuuL-c- 
-9 .9590 504 

U -2 0308503 
4.28576-03 
,5.0006503 



PR NS, linear model, Rop=0.34 
Re= 15,000, at the 16,000 time step 
NCC-1 1 4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 0, ctt=4 Oe-4 
Convergence 3 order or Max iter=120 



■ 


2. 90 18 &02 
2.7051502 
2 EO07F-Q2 
2.3122602 
2 11386-02 
l 9193L-02 
1.7229502 
1.52W502 
1 3300FQ2 
1 1335502 
9.3700503 
7 .4002503 
54417503 
3.4772503 


PR NS, linear model, Rep =0.34 
Re= 15,000, at the 16,000 time step 
NCC-1 1 4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max iter=120 



I 

i 


mu 

3.3901504 
3 1613504 
2 9326F-04 
2.70306-04 
2 4731604 
2.24936-04 
2.0170504 
1 70005 04 
1 SOOIf EM 
1.3314504 
1 1O205EM 
0 7309505 
6 4515505 
4.1640505 
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B.5 PRNS with RCP = 0.4 


History of velocity components at probe 1 and probe 14: 




Pipe flew Re=15k, at probe 1 4 
PRNS linear model, RCP =0 38 
2nd=Q.Q. 4th=Q.Q0GQ1 , cfi=1.0, dt=4.Qe-4 
Convergenceorder=3. Max iter=120 

u 


Pipe flew Re=15k, at probe 1 4 
PRNS linear model, RCP =0 38 
2nd=Q.Q. 4th=O.QOOQ1 , cfl=1.0, dt=4.Qe-4 
Convergence order=3. Max rter=120 
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History of subscale turbulent kinetic energy at probes 1 and 14. 




Contours of w, u, v, pg, k, ju + p T at the time step 20,000: 


PRNS, linear model, Rop-0.3 Q 
Re=15000, at the 20,000 time step 
NCO-1 1 4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=Q .00001 , CFL=1 0, ctt=4 Oe-4 
Convergence 3 order or Max tter=120 



i 


1 .2444&-Q0 

i.igsoe+oo 

1 n73fiir*nn 
y utnoL-o! 
9 O3T0E-P1 
U 1737601 
7.31076O1 
S.40506Q1 
SfilKSFdl 
4.TST5&01 

3.9035601 

3.043411-31 

2 1954601 
1 3413&01 


PRNS, linear model, Ftap=0.3£ 

Re=15000, at the 20,000 time step 
MCC-1 .1 A. no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



1 5954605 
1. 3544605 
1 l?3.ir Cf. 
0.9223609 


5 511*60(5 
4.3010600 


■ -3 20345,07 
*' -3 (33 ID I? EK5 
-*.SH17L-0(i 
-7 2523600 

1 - 9.5930600 
.1 10746.05 
■ 1.4194605 




PR US, linear model, Rop»Q,36 
Re=1500Q, art the 20,000 time step 
NCO-1 .1 A. no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 O, ctt=4 Oe-4 
Convergence 3 order or Max tter=120 



| 4.3951602 
,71 99602 
1 3W37F-Q3 
' 2.399564)2 
1 9904602 
1.0142602 
3 .300 1603 
I 4.3313603 
1 0 143603 
-1 (jyUtiL-02 
-2.3007602 


I 


4.0429602 
4.7190602 
■4.39526 02 


PR NS, linear model, Rep =0.38 
Re= 15000, art the 20,000 time step 
NCO-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, cit=4 Oe-4 
Convergence 3 order or Max tter=120 



P9 

1.5050602 
1.3470 602 
1 1090602 
1 U3 17602 
h T 392603 
7 1554603 
5.5747603 
3,0930603 

?4i3?raa 


0.3340604 

-7.4030604 


1 -2.3291603 
-3.0009603 
-5.4000603 
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PR MS, linear model. Rep =0.33 
Re= 15000, art the 20,000 time step 
NCO-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



I 


2.0010E.O2 

2.0121&O2 

2.23:213 £*02 


■ CW29E-07 

1 0532E-O2 
1 003S&O2 
1 4707&O2 
i :Mnrn7 
I UW3L-02 
3.CH50E-03 
7.U04&O3 
5.201 1&03 
3 3530&Q3 


PR MS, linear model. Rep =0.33 
Re=1500Q, art the 20,000 time step 
MCC-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, ctt=4 Oe-4 
Convergence 3 order or Max tter=120 



I 

i 


mu 

3.5007&O4 
9.35S0&W 
3 1 1 14 F- EM 

2 0377&O4 
7 m i&w 
2 .3005 £*04 
2.1300&O4 
1.0832&W 
1 PJ40SE EM 
1 JQ53E-EM 
1 I022E-O4 

3 1000605 
0.7490&OS 
4.3132&05 



B.6 PRNS with RCP = 0.5 

History of velocity components at probe 1 and probe 14: 



Pipe flew Re=15k, at probe 1 
PRNS linear model, RCP=0,5O 
2nd =0.0. 4th=0.Q00Q1 , cfl=1.0, dt=4.0e-4 
Converged ceorder=3. Max iter=120 
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Pipe flow Re=15k, at probe 1 4 
PRNS linear model, RCP=Q 50 
2nd =0.0. 4th=0.QQ0Q1 , cfl=1.0, dt=4.0e-4 
Convergence order=3, Max iter=120 

u 


Pipe flow Re=15k, at probe 14 
PRNS linear model, RCP=0 50 
2nd =0.0. 4th=0.00001 , cfl=1.0, dt=4.0e-4 
Convergence order=3, Max iter=120 



- 0.15 - 


I ■ ■ ■ ■ I ■ ■ ' ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I - i - - 1 ■ i ■ ■ I 

0 2500 5000 7500 10000 12500 15000 

Time steps 


i ■ ■ ■ ■ i ■ i ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i - i - - 1 i i ■ ■ i 
0 2500 5000 7500 10000 12500 15000 

Time steps 


History of subscale turbulent kinetic energy at probes 1 and 14: 




Contours of w, u, v, pg, k, ju + p T at the time step 15,000: 


PRNS, linear model. Rep =0.50 
Re= 15000, at the 20,000 time step 
NCO-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max iter=120 



w 

1.2514&-Q0 

1.1059*+DD 

i nana^on 
9.9400*01 
a [1920*111 
0 2375*01 
7. 3023*01 
*5270*01 
S 6710P (11 


-J 0105*01 
3.O013E-C1 


I 


3 IU0UL-01 
2.2500*01 
1.3955*01 


PRNS, linear model. Rep =0.50 
Re= 15000, at the 20,000 time step 
NCO-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=Q .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



1.2 ISO 60S 
1 .0179*05 

aragar-Qa 

0.2397*08 
4 2003*08 

2 2300*00 

3 2707&O? 
1 *127500 


•3 013nF00 
-5.5633*00 
-7.5537*00 
-3.5240*00 
-1 1434*05 
• 1 9405*05 
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PR NS, linear model. Rep =0.50 
Re= 15000, art the 20,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0.0, 4th=0 .00001 , CFL=1 0, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



3 1Q1T&02 
2.02456O2 
2 1472602 


1.0700602 
1 1827602 
7 1540603 
2.3023603 
■2 .3903603 
.7 1920603 


I 


-1 1935602 
-1 .0709602 
-2 140U6EI2 
-2.0253602 
-9.1025602 


PR NS, linear model. Rep =0.50 
Re=15000, art the 20,000 time step 
NCC-1 .1 .4, no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



I 


1.5570602 
1.4009602 
1 2439 F- 02 
1 .00706 02 
9 3000603 
7.7313603 
0.1020603 
4 5927603 
3 0234F 03 
1 4541603 
-1.1519604 
-1 004S6O3 
-3.2530603 
4.9231603 



PR NS, linear model. Rep =0.50 
Re= 15000, art the 20,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce= 1 00290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



I 


1U13V 

2 7509602 
2.5993602 
2 3010602 
2 1340 602 
2 0OM&O2 
1.3100602 
1 0312602 
1 .4435602 
1 2559 r 02 
1 0003602 
0 0009603 
0.3307603 
5.0545603 
3.1793603 


PR NS, linear model. Rep =0.50 
Re= 15000, art the 20,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce= 1 0D290 

2nd =0 0, 4th=0 .00001 , CFL=1 O, dt=4 Oe-4 
Convergence 3 order or Max tter=120 



4 4099604 
4.1743604 
1 3 0S70F-04 
3.5510604 
3 2543604 
2.0470604 
2 8410604 
2.3343604 
2 027(5 F 04 


1.7210604 
1 4143604 
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1 1077604 
0.0 100 605 
4.9434605 
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B.7 PRNS with RCP = 1.0 


History of velocity components at probe 1 and probe 14: 




Pipe flew Re=15k, at probe 14 
PRNS linear model, RCP- 1.0 
2nd=Q.Q. 4th=G.QQ0Q1 , cfl=1.0, dt=4.0e-4 
Convergence ordef*3, Max rter=120 

u 


Pipe flew Re=15k, at probe 14 

PRNS linear model, RCP=1.D 

2nd =0.0. 4th=G.Q00Q1 , cfl=1.0, dt=4.0e-4 

Convergence order=3, Max rter=120 
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History of subscale turbulent kinetic energy at probes 1 and 14: 


Pipe flew Re=15k, at probe 1 
PRNS linear model, RCP=1 0 
2nd=Q.Q. 4th=G.QGQQ1 , cfl=1.0, dt=4.Qe-4 
Convergence order=3. Max iier=120 

k 


0J03 - 


0.02 - 



-am - 


• I 1 ■ ■ ■ i I 

0 2500 5000 7500 10000 12500 15000 

Time steps 



NAS A/TM— 2008-2 15418 


73 


Contours of w, u, v, pg, k, ju + p T at the time step 15,000. 


PRNS, linear model, Rcp-1.0 
Re-15000, at the 15,000 time step 
NCO-1 1 4. no smoothing 
P-referen ce= 1 00290 

2nd -0 0, 4th-0 .00001 , CFL=1 0, dt-4 Qe-4 
Convergence 3 order or Max tter=120 



I 


1 12Q7E+0D 
1 .04456*00 
9 002B IT- 01 
0.9207601 
0 1S00E-D1 
7 .3993601 
0.0344601 
5.0723601 
S 11036-01 
4 .3402601 
3.5001 6-01 
2.0240601 
2.Q010&O1 
1.2999601 


PRNS, linear model, Rcp-1.0 
Re-15000, at the 15,000 time step 
NCO-1 .1 A. no smoothing 
P-referen ce- 1 00290 

2nd -0 0, 4th-0 .00001 , CFL-1 0, dt-4 Qe-4 
Convergence 3 order or Max £ter=120 
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Convergence 3 order or Max 6ter=120 
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PRNS, linear model, Rcp-1.0 
Re-15000, at the 15,000 time step 
NCO-1 1 4. no smoothing 
P-referen ce- 1 00290 

2nd -0 0, 4th-0 .00001 , CFL-1 0, dt-4 Oe-4 
Convergence 3 order or Max iter=120 
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PRNS, linear model, Rcp-1.0 
Re- 15000, at the 15,000 time step 
NCO-1 .1 .4, no smoothing 
P-referen ce- 1 00290 

2nd -0 0, 4th-0 .00001 , CFL-1 0, dt-4 Oe-4 
Convergence 3 order or Max tter=120 
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PRNS, linear model, Rcp-1.0 
Re- 15000, at the 15,000 time step 
NCC-1 .1 .4. no smoothing 
P-referen ce- 1 00290 

2nd -0 0, 4th-0 .00001 , CFL-1 0, dt-4 Oe-4 
Convergence 3 order or Max tter=120 
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Appendix C: Effects of anisotropy and rotation in the subscale model 


C.l Nonlinear PRNS with RCP = 0.3 at Reynolds number 15,000 


History of velocity components at probe 1 and probe 14: 



Pipe flew Re= 15 k, at probe 1 
PRNS non-linear model, RCP=D .30 
2nd=Q.Q. 4th=G.QG1 , cfl= 1 .0, dN4.0e-4 
Convergence order=3, Max iter=120 
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History of subscale turbulent kinetic energy at probes 1 and 14. 


Pipe flaw Re=15k, at probe 1 
PRNS non-linear model, RCP=D.30 
2nd=Q.Q. 4th=0.Q01 , cfl= 1 .0, d\=4.Qe-4 
Convergence order=3. Max rter=120 


k 


0.002 - 




0.001 - 


o - L 

— L 

0 




5000 7500 10000 12500 15000 



Time steps 


Time steps 


Contours of w, u, v, pg, k, ju + p T at the time step 15,000: 


PRNS, non-linear model, Rep =0.30 

Re=15,000, at the 15,000 time step 

NCC-1 1 4. no smoothing 

P-referen ce= 1 0D290 

2nd=0 0, 4th=0.Q01 . CFL=1 dt=4.0e-4 

Convergence 3 order or Max iter=120 
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PRNS, non-linear model, Rep =0.30 

Re= 15,000, at the 15,000 time step 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 0D290 

2nd =0.0, 4th=Q.Q01 . CFL=1 .0, dt=4.0e-4 

Convergence 3 order or Max iter=120 
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PRNS, non-linear model, Rep =0.30 

Re= 15,000, at the 15,000 time step 

NCC-1 .1 A. no smoothing 

P-referen ce= 1 0D29G 

2nd =0 0, 4th=0 001 , CFL=1 .0, dt=4 Oe-4 

Convergence 3 order or Max iter=120 
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PRNS, non-linear model, Rep =0.30 

Re= 15,000, at the 15,000 time step 

NCC-1 .1 .4, no smoothing 

P-referen ce= 1 0D29G 

2nd =0 0, 4th=0 001 . CFL=1 .0, dt=4 De-4 

Convergence 3 order or Max iter=120 
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PR NS, non-linear model, Rep =0.30 

Re= 15,000, at the 15,000 time step 

NCC-1 .1 .4. no smoothing 

P-referen ce= 1 00290 

2nd =0.0, 4th=0.Q01 . CFL=1 .0, dt=4.0e-4 

Convergence 3 order or Max iter=120 
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PR NS, non-linear model, Rep =0.30 

Re= 15,000, at the 15,000 time step 

NCC-1 .1 .4. no smoothing 

P-referen ce= 1 0D290 

2nd =0 0, 4th=0 001 , CFL=1 .0, dt=4 De-4 

Convergence 3 order or Max tter=120 
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C.2 Nonlinear PRNS with RCP = 0.3 at Reynolds Number 150,000 
History of velocity components at probe 1 and probe 14: 


Pipe flow Re=1 50k, at probe 1 
PRNS, RCP=0.30 

2nd =0.0. 4th =0.001 , cfl=1.0. dt =2.5^5 
Convergence order=3. Max iter=120 
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Time steps 


Pipe flow Re=1 50k, at probe 1 
PRNS, RCP=0.30 

2nd =0.0. 4th =0.001 , cfl=1.0, dt =2.5^5 
Convergence order=3. Max iter=120 
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Pipe flow Re=1 50k, at probe 1 4 
PRN5, RCP=0.30 

8 r 2nd =0. 0. 4th =0 .00 1 , cfl = 1 .0 . dt -2.5^5 
C on vengen ce order=3 . Max iter= 1 2Q 
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History of subscale turbulent kinetic energy at probes 1 and 14: 


D.5r Pipe flow Re=150k, at probe 1 
PRNS, RCP=Q.3Q 
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Convergence order=3 r Max iterM20 
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Contours of w, u, v, pg, k, /u + ju T at the time step 38,400. 


PR NS, non-linear model, Rep =0.3 

Re=1 50000 , 36400 time steps from linear model 

NCO-1 .1 .4, no smoothing 

P-referen ce= 1 00290 

2nc=D.O, 4th=0.001, CFL=1.0, dt=2 5e-5 

Convergence 3 order or Max tter=120 
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NCO-1 .1 4. no smoothing 
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Convergence 3 order or Max iter=120 
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